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AoCTRACT 


Tko. 6tfLactaA.t Oy a nofiraaZ ^kodz v>avt Zn. rnon-atoraZc coJiQ.6 cd 
liZQk i'ladt numb^fi 6 vjhiin Zon.Zzo.tZon. '"“Cconc-A ZmpoKtant -Li tk^ iabjadt 
£?{,' t.LZ p^zitnt itudj. razi iM 30 iZnce. tht v)Zdth o{s tkt fitZ.axcitZon. 
zone Za cor,v..azia.bZ^ to .tht yoi djncriZa ihodz r.jZdtk a. iZp.attane.oa 6 
itadj o-Q tilt tvjo KzcZoYLi Zi> n.aqaZfic.d. TZziour'i ri izZmtZc .tkzoa:j 
avpzLoack uiZn^ t'lz Zottziiann e.C[acitZon c it-i oi e.qaatZ.OYii n.tpn.z- 
ie,ntZn'j the entZ^e non-eqaZtZb.xZarn ztcgZon have been deaZved 
i/'>hZch atie iotved nap.eKZeo.tt. 3 . nadZatZon and mattZpte ZonZzo.tZon 
eileeti have been Zgnozted. 

A tziZmodat dZ 6 t-xZbatZon ^anetZon lioi been aied to A.epaei ent 
the atom gai and M.axuJettZan dZitziZbutZon ^anetZom {^oa Zon and eteet- 
gaiei. i^taitZc cottZiZoni betvjeen the dZ'^-^zaent ipeeZei have been 
eoYiiZdefied tkJtoagk Znveoi e-poojen. tavi io/iee mod.eti. cxpeaZmentattij 
deteamZned ZnetaitZc eottZiZon eaoii ieetZons hove been aied {^on. 
atom-atom and atom- eteetzion ZnetaitZe cottZiZom. The coaaeipond- 
Zng aecombZnatZon aatei have been outaZned ij-xom the Coha. equatZon. 
Fax.thex.mox.ej an addZtZonat -moment 0 ^ the Cottzmann eoaatZcn hc..i 

X 

been dex-Zved.. 

UiZng the appx.opx.Zo.te hoandox.'j eondttZom the x.eiattZng eqaa- 
tZoni together ujZth the CoJio eqaatZon gZve the jump eondZtZom 
aaxoii the ihock uave and have been evc.taated namex.Zeotltj , A name- 
xZcat iotatZon oi the ihock itx.actax.e eqaatZoni on the ICM 7044 
compatex hai been obto.Zned (iOX the ihock wave pxoiZtei ^ox axgon, 
kxypton and xenon gaiei i^ox a xange o{. Mach namoexi ixom 15 to 40 



and ^oA a Aangt LnLt-lat toM-po^AatixAdh and pAd^&uAz^ , In addZ- 
t-ion e.ZdctA'icai. c.ondu.atAv^ty and AC'I^Aact^vd -index -in the -ihock- 
o’javd have been aaZcuZated, The -ikock th-icknd.&.6 and the tk-ickm-6.6 
paAametdA tu dti-ined. b>j fe-t-bckek and. utjAon hove been computed 
and a. good agAeernent o{- the tktdznehi, panametea t& lound vitth 
the expeatmentaZ Ae-6uZtA o£ '‘long and ocAi-hadeA. 
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CHAPTER 1 


IMTROPUCTIOIT 


1 . 1 G-e^e raH j 

The structure of a shock wave is of interest for 
various reasons. In the beginning it attracted the attention 
of many scientists as a theoretical problem, the solution of 
which aided the understanding of the physical mechanism of 
shock compression, a remarkable phenomenon xn gas dynamics. 

At very high altitudes where the densitjr of air is very low, 
the mean free path may become comparable in size to the body 
dimensions. As the shock thiclcneos is only of the order of a 
few mean fiee paths the variation of the macroscopic quantities 
like pressure, temperature, velocity etc. in this region are 
of great interest. When it was realised that the hydrodynamic 
equations of motion lose their significance in representing 
the phenomenon inside strong shock waves which are of the order 
of a few rneai free paths, attention was diverted towards the 
kinetic theory approach to the fluid dyna''iical problems. The 
one-dimensional shock structure problem duo to its simplicity 
and freedom from boundary efiects attra.cted the attention of 
scientists working in kinetic theoiy of gases as a very good 
example for checking the validit 3 r of the approximate solutions 
proposed for the Boltzmann equation, later shock waves vrere 
produced in laboratories to reaHise high temperatures aid study 
the various high temperature phenomena occurring in gases like 
excitation of various degrees of freedom, dissociation. 
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chemical reactions, ionization etc. A comparison of 
experimental results with theory facilitated deducing much 
Yaluahle information about the reaction rates of these 
processes. The study of shock waves of veiy high strength 
in which radiation plays an important role helps in clarifying 
the various interesting asi'iccts like shock front luminosity 
and other optical effects oId served in strong esrplosions in 
air. Shock waves of very high Mach numbers are also of 
interest in electro -magnetic shock tube work. With its 
far-from-equilihrium properties, its sim£jle geometry and 
boundary conditions, the shock wave scructure prohlem is a 
unique choice for testing non-equilihriuii theories. 

1.2 She c_k J(ave^s__ in a Wpji re act in^j^ 

The equations of conservation of mass, moiiiGntum and 
energy for a one-dimensional flow of an ideal gas admit a 
non-trivial solution which represents a possible finite 
jump in the flow quantities like velocity, temperature, 
pressure, in the absence of a characteristic length scale for 
the problem. Associated with this jump in pressure, velocity 
and temperature is an increase in entropy which indicates 
a dissipation of mechanical energy and an irreversible 
conversion of mechanical energy into heat. In order to 
understand this transition it is therefore necessary to 
include the dissipative mechanisms viz. viscosity and heat 
conduction attributable to the molecular structure of the fluid. 
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These processes gi’ve rise to additional transfer of momentum 
and energy resulting in nonadiabatic flow and irreversible 
transfer of mechanical energy into heat. With these dissipative 
mechanisms the discontinuities in the flow quantities appear 
as steep gradients. The shock wave, a plane of discontinuity 
for ideal fluids, becomes a wave of finite thickness. Due to 
the large gradients in the flo?/ parameters existing in this 
region, the dissipative mechanisms are important only inside 
this region whereas outside the shock they are ineffective. 

Even though the flow pa-rameters change scrictly over the range 
-oo<x<+‘=<= a large change occurs only in a very small region of 
the order of a few mean free paths. Equilibriuiv can be assumed 
to exist outside this region of large gradients. Though the 
transport properties, viscosity and thermal diffusivity are 
similar in nature their relative roles in the shock wave are 
not alike. Viscosity plays the role of converting the directed 
molecular motion to random motion thus converting kinetic 
energy of mass motion to heat energy inside the shock wave. 
Conduction plays the role of redistributing the energy of 
random motion or heat in space and hence in heat transfer from 
region of higher temperature to lower temperature. The entropy 
change associated with viscosity is due to the generation of 
heat from mechanical energy and hence is positive everywhere. 

On the other hand the change in ent.ropy associated with heat 
conduction is due to net inflow or outflow of heat at any point 
which can be positive or negative. Thus associated with heat 
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conduction a local decrease of entropy is possible. Howe-ver, 
the overall entropy change across shock due to heat conduction 
is positive. 

The solution of the one-dimcnsional hydrodynamic 
equations neglecting the viscosity hut retaining heat conduction 
studied hy Raleigh, revealed that a continuous solution for the 
flow quantities is possible for weak shocks whereas such a 
continuous solution may not be possible in general for strong 
shocks. Indeed, it is impossible to construct a continuoas 
solution for the flow parameters far a sufficiently strong shock 
wave. On the other hand it is possible to construct c.ntinuous 
solutions of all flow parameters for all shock strengths by 
considering viscosity even if the heat conduction is neglected. 
This difficulty, explains the importance of viscosity in achieving 
the irreversible compression of the fluid and that a strong 
discontinuity can disappea,r only through the action of viscosity 
but not of heat conduction. In this 'respect effect of 
conduction is secondary in the shock wa.ve. Studies involving 
only thermal conduction are of interest since they illustrate 
the effect of other heat transfer mechanisms like radiation, 
electron heat conduction in plasma, etc. Thus in studying the 
flow variables inside the shock it is necessary to include the 
effect of viscosity, primarily responsible for smoothening the 
discontinuity in flow quantities. 
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1 . 3 Th e Probl ein p^f _Sh q ck^ _S t m c t ur e : 

The fact that the flow is alaost in equilihriun both 
upstream and doTOstream of the shock and that gradients in the 
flow quantities are negligihle in these regions implies that 
the flow parameters on the downstream side can he related to 
those of upstream side hy simple hydrodynamic equations known 
as Rankine-IIugoniot equations for the shock front. These 
relations give the correct values of the downstream parameters 
as long as there is no change in the microscopic structure of 
the gas as happens at high temperatures. The problem of shock 
structure y;hich does not involve a change of microscopic 
structure of the gas (excitation of internal degrees of freedom 
and a breakup of molecules) boils down to a study of the 
transition region in which the viscosity and thermal conductivity 
are active, boimded on either side by regions of equilibritim 
whose parameters are related by R-H relations. The fact that the 
flow parameters on the downstream side do not depend upon the 
manner in which they vary in the transition region demands that 
any theory ' proposed to explain the mechanism inside the 
transition region should give the downstream values as 
predicted by the R-H relations. When the internal degrees of 
freedom are excited those Rankine-Hugoniot relations have to 
be suitably modified to yield appropriate downstream equilibrium 
flow parameters. 
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"1 • 4 - Stock Wave Tlaickn_e_ss : 

A simple and suitable parameter is to be identified 
for comparing the various theories and experiment s. An obvious 
choice of this parameter is the thickness of the transition 
region (if one such can be defined) due to its simplicity. The 
definition of such a thiclmess is inva>,ri0,bly arbitrary if one 
wants a simple, relevant, representative i-janameter. Since the 
density, flow velocity, pressure and temperature do not vary in 
the same manner each can be Lised to define a thickness. These 
thicknesses need not be equal though they may be of the same 
order of magnit-ude. The arbitrariness also arises from the 
definition of^a thickness’ for any of these profiles, because 
these quantities vary continuously in <x<+°o. However, as 
the major variations of the flow quantities are confined to a 
narrow region, a suitable definition of s^^ch a thicloiess can 
not be too far from realityi one such definition is, the region 
in which a flow variable changes by an ’accepted fraction’. 

But the thickness is veiy sensitive to the choice of this 
’'accepted fraction.' 

The most widely vised definition, however avoids this 
ambiguity and replaces the actual profile by a linea,r profile 
with a slope equal to the rao.ximum slope of the actual profile. 
The thickness is that region in which any chosen flow parameter 
varies linearly from the upstream value to the downstream 
equilibrium value. The major objection to this definition 
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is that it is too localized in description and fails to make 
sense when the profile han a point of inflexion with T/ertical 
slope. This definition is not suitable 'when the slopes of the 
actual profile at distances away from mazimum slope point are 
■very much different from ma,xtmum slope as v/hen ionization 
occurs ‘behim’. the shack. xTewertheless, for the description 
of conventional shock v/aves ass-uming that the flow parameters 
yaryc smoothly, the above definition is suitable since it can be 
evaluated very often without undue difficulty. Other 
definitions such as Grad's area rule are also often chosen. 

It would seem reasonable to choose thar definition which agrees 
well with the experimental mea.surement s; unfo’rtunately the 
experimental values themselves involve a. certain amount of 
uncertainty depending upon the type of the experimental 
technique . 

lor shock waves involving ionizalion, Petschek and 

62 77 

Byron and also Wong and Bershader defined a thickness '6' 

as the distance behind the gasdynamic shock front to the point 

where the degree of ionization reaches 70^ of its downstream 

equilibrium value. Por purposes of comparison, they defined 

a parameter equal to the product of the upstream equilibrium 

pressure and the elapsed time in the laboratory reference frame 

between the passage of the shock front and the 70‘% ionization 

point . 
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* 5 Lengt h f^r th_e_^ _Stip_ck J^e : 

It is important to stress the necessity of a length 
unit in the problem of shock wa've structure to nondimensionalise 
any macroscopic length scale of the problem. The absence of a 
body (which would have provided a characteristic length) in the 
problem poses considerable difficulty in identifying one such 
imit of length. The kinetic theory description of the gas flow 
suggests the 'mean free path’ a,s a suitable parameter, suitable 
because the molecular collisions are responsible for the 
relaxation or equilibration. The mean free path has often 
been used as a length scale for the shock structure problem even 
when Navier-Stokes equations describe the motion. In the latter 
case mean free path is calculated based on the measured values 
of some macroscopic property of the gas described through the 
kinetics of molecular collisions. One such ’mean free path' 
used quite often is the viscosity mean free path, for a given 
viscosity a mean free path can be defined for any gas so that 
measurement of the viscosity is equivalent to the measurement 
of the viscosity mean free path for the gas. Thus in the 
continuum theoi^y there is no difficulty in defining the mean 
free path, but the problem is wdiere that mean free path should 
be measured. 

The problem is a little more involved in the kinetic 

theory approach since the mean free path depends on the inter- 
molecular force law chosen, further very often a criterion 

is to be chosen to define a ’collision' or a collision cross 
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section since in the strict sense the col3.ision cross section 
is infinite in such cases. Quite often the -viscosity cross 
section is chosen for this purpose. Further since it is not 
al-ways possible to perform the calculations with the true force 
lawsjmethods have been devised to use the knovoi viscosity- 
temperature dependence so that valid comparisons can be made 
with real gases. Even so ’where’ the mean free path is to be 
evaluated still remains open. 

In the case of the relaxation region where the effect 
of inelastic collisions is dominant and viscosity plays a 
negligible role, it is appropriate to use the mean free path 
CO nre spending to the inelastic collisions. In cases when there 
are more than one type of inelastic collision, as in the case 
of ionization due to atom-electron and atom-atom collitions, 
recombination collisions etc. it may be proper to chaocc a 
mean free path corresponding to the dominant collision. 

As regards ’where ’the mean free path is to be measured, 

the upstream mean free path has been used -very often as it is 

simple to calculate. The choice of the upstream mean free path 

makes the nondimensional shock thickness depend strongly on the 

choice of the inter-molecular potential at large Mach numbers, 

55 

Narasimha et al suggested a mean free path corresponding to 
the hot side as a proper choice which suppresses the sensitivity 
of the mean free path to the inter-molecular potential. The 
upstream mass velocity can also be advaiitageously used to 
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construct a length scale in place of the thermal speed. Such 

58 

a length has been used hy Oherai , even though the advantage 
of it has perhaps not been rea.lised. This length scale has been 
used in the present investigation to nondimensionalise the 
equations for the entire region. 

^ The Relaxa tion Phenomena jL_n jla.se 

In the case of veiy strong shock waves in a real ga.s, 
the high temperatures a,ttained result in the excitation of 
various degrees of freedom. The state of such a gas depends on 
the concentration of the various components such as molecules, 
atoms, ions and electrons, and also on the distribution of the 
internal energy among the various degrees of freedom. Generally 
the internal energy of a gas consists of the energy of 
translational motion of the particles, rotational and vibrational 
energies, chemical energy of reactions, ionization energy, and 
electronic excitation energy of the molecules, atoms ajid ions. 
Associated v/ith the excitation and reestablishment of 
thermodynamic equilibriiom of any of these degrees of freedom 
there are corresponding chn,racteristic times known as the 
relaxation times. Relaxation times for exciting the various 
degrees of freedom differ appreciably from each other so that 
under ceirtain conditions it may be possible to achieve 
thermodynamic equilibrium for some but not all degrees of 
freedom. Since the kinetic energy and the momentum transferred 
during a collision between particles of comparable masses are 
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on the average of the same order as the kinetic energy and 
momentum of the colliding pa-rticles, the relaxation time for 
establishing ecLuilihrium of translattionaJ degrees is of the order 
of the average time between gas kinetit; collisions. Since the 
collision times are small compared to the flow times of 
manposcopic interest, it is possible to ascirlbe to the ga,s at 
every instant of time a 'translational' temperature, which 
characterizes the average kinetic energy of translational motion 
of the particles. Under conditions of partial thermodynamic 
equilibriumCthermodynamic eq.uilibrium for the individual degrees 
of freedom) it is implied that the distribution of energy and 
concentration of the respective components of the mixture for 
these degrees of freedom is in equilibrium with the translational 
temperature of the gas. But quantities associated with 
nonequilibriuim degrees of freedom may be arbitrary depending 
on many factors including the history of the gas. Such is the 
situation encountered in a shock front. The distribution of 
energy concentrations is not determined in this case simply 

by the temperature, density and composition of the gas as in 
the case of thermodynamic equilibrium, but also by the kinetics 
of the physico-chemical processes leading to the reestablishment 
of equilibrium for the given degrees of freedom, 

formally the establishment of equilibrium of various 
degrees of freedom follows the following hierarchy translational 
degrees, rotational degrees, vibrationa.! degrees, dissociation 
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and chemical reactions, ionisation and electronic excitation. 

Due to the pronounced difference in the va.rious rela:^ation times, 
it is vexy often possible to study some relaxation processes 
separately, isolating them from the remaining ones. 

The study of the kinetics of physico-chemical 
relaxation processes has two parts: the first dealing with the 
rates of the elementary processes leading to the excitation 
of one or other degree of freedom, the second dealing with the 
kinetics of the rela.xa,tton process itself under specific 
conditions, sls far as the rota,tional degrees are concerned, 
it is practically alv/ays possible to assume that equilibrium is 
established as quickly as with the translational degrees of 
freedom so that they can be considered as classical type. On 
the other hand the problem of vibrational relaxation becomes of 
practical importance when the vibrational energy of the molecule 
is comparable to the translational energy in which case the 
vibration mode takes an essentially quantum character. When 
the translational energy is small companed to the vibrational 
energy of the molecule the excitation of the vibrational 
degree can be ignored whereas in the other limit of far 
classical region (the translational energy is large compared 
to the vibrational energy of the molecule) almost all the 
molecules will have dissociated and the pinblem loses its 
meaning. The dissociated gas behaves as a monatomic gas in 
which only ionization and electronic excitation are possible. 



13 


1*7 izi3a£_ Sho^ck_^ Mo^ gmi c, 

In studying the ionization and its effect on shock 
structure it is interesting to deal with a monatomic gas. Due 
to the absence of a large number of degrees of freedom possessed 
by polyatomic gases it is easier to obtain higher temperatures 
in monatomic gases. Monatomic gases are also simpler since 
apart from electron state excitation ionization is tho only 
relax;ation process responsible for broadening the shock front. 
Most of the studies have been made with noble gases argon, 
krypton and xenon because their low ionization potential results 
in a higher degree of ionization for a given shock strength. 

Also their greater atomic weight (lovrer speed of sound) helps 
the production of strong shock waves. Argon has been used most 
often because its atomic propeities a.rG v/ell known and the pure 
gas is inexpensive, 

lo study the physics of an ionizing shock wave, a 
stationar7 strong shock wave in a monatomic gas can be considered 
for simplicity. Upstream of the shock the neutral gas is in 
theimodynamic equilibrium. As the particles travel across the 
shock their ordered motion will be disturbed and a large part of 
the kinetic energy of the flow will be converted into thermal 
energy of random motion as in an inert gas shock. As a result 
of the high temperatures developed there will be a large number 

of collisions in which energies of the order i (ionization' 
potential of the atom) will be transferred to the electrons in 
their orbits resulting in ionization of the atom. The 
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electron- atom collisional ionization is -very effective due to 
its large cross section. However, this process needs some 
initial priming electrons to form an electron avalanche. One 
of the mechanisms which can lead to this initial ionization is 
throtigh atom-atom collisions. As the cross section for this 
atom-atom collisional ionization process is extremely small it 
takes a large time for the formation of priming electrons. 
Conseq.uently 5 significant ionization takes place only downstream 
of the gas dynamic shock with its structure determined hy the 
elastic collisionaA dissipative mechanisms as in an inert gas. 
However, at Mach numbers greater than -^20, vifhen the temperature 
inside the gas dynamic shock is high enough to cause significant 
ionization by atom-atom collisions, the structure of the gas 
dynaflic shock front is appreciably influenced by the ionization 
process. Hence inelastic collisional processes should be 
included in the study of the gas djnax-iic shock also. Other 
possible mechanisms of producing priming electrons are 
photoionization atid ioniza-tion due to impurities. 

Avalanche ionization begins when the rate of 
atom-electron collisional ionization becomes greater than the 
rate of the initiating process. Since the rate of atom-atom 
collisional ionization is very low, avalanche ionization begins 
with only a very few electrons. The point at which atom-electron 


ionization rate overtakes the atom-atom ionization rate is 
referred to as the 'take over point' or the 'onset point'. 
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The region dowistrean of the onset point is relatively 
insensitive to the atom-atom ionization cross section. At a 
constant electron temperature 5 the avalonche increases 
exponentially with time (disTance) until recombination begins to 
balance the ionization rate lea,ding to equilibrium. In actuality 
the avalanche proceeds in a much more complex manner. 

Pirst, the ionizatiqn occurs through a two-step process in which 
an atom gets excited to its first excited level and then 
subsequently ionized, further ea.ch ionization results in the 
electron gas losing an- amount of energy I, corresponding to 
the thermal ^ergy of many electrons so that the electron gas 
tempera,ture is lowered. If this energy is not replenished 
the electron gas temperature would continue to fall veiy 
rapidly, A fa,ll in the electron temperature in turn lowers the 
ionization rate which depends very strongly on the electron 
temperature through Boltzmann factor e“^/^e . The electron 

energy loss resulting from ionization will be compensated by 
energy transfer from the hightenp^'roturu heavy particle gas 
through ol.asti© collisions. However the energy transfer 
between the heavy particles and electrons proceeds very slowly 
due to the large difference in their masses restricting the rate 
of electron-atom collisional ionization to a great extent. The 
energy an electron gains during a collision with an ion is 
several times more than that with an atom due to the coulombic 
nat-ure of the force field between charged particles, nevertheless 
in the initial stages when the ion population is ve^y small. 
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the at ora- electron elastic collisiona-1 energy transfer is still 
significant. Because of the nearly equa,l masses of atoms and 
ions and also the large elastic collision cross sections 
Between them, energy exchange Between then proceeds very 
rapidly. Consequently the atoms and ions have the same . 
tempera>,ture throughout. Thus^But for an initiad. region, ions 

r 

very efficiently supply energy to the electron gas serving as 
an intermedia,ry in the energy flow from the atoms. In the 
electron gas the energy is very ra.pidly distriButed through 
elastic collisions so that an electron temperature can Be 
assigned which is differeiit from the heavy particle gas 
temperature. In the electron- atom collisiona.1 ionization 
zone energy lost By the electron gas is Balanced By the energy 
gained through elastic collisions with heavy particles so that 
the electron gas temperature changes very little. 

Further downstream when the electron population is 
Built up close to its equiliBrium value the rocoiiBination in 
which ions and electrons unite to forr-i neutral a,toms Becomes 
important. This recoriBination rate Balances appreciably the 
ionization rate so tha,t the net rate of electron production 
decreases. In this region the atom temperature starts falling 
finally approaching the electron temperature when thermal 
equilibrium is reanhed. later, the electron density approaches 
its equilibrium value at which the recombination rate exactly 
balances the ionization rate and chemical equilibrium is 
reached. 
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® P olarization in sn Ionize _d G-a.s : 

If an electron density gradient exists in a gas, the 
electrons tend to diffuse so as to smooth out the gradient, due 
to the high mobility of the electrons resulting from their 
exceedingly small mass. The ions, hov/over, have a much lower 
diffusion velocity due to their ls,rge mass, and, therefore, 
cannot move relative to the neutral atoms. This leads to the 
motion of electrons relative to the ions resulting in a charge 
separation and creation of a strong electric field with its 
intensity given by Gauss equation. This electric field 
however prevents further polarization and inhibits the 
diffusion' electron current. 

Since the field due to a charge in a plasma is 
shielded within one Debye length the electrons cannot be 
separated from ions by more than this distance. If the Debye 
length is much smaller than the electron mean free path, 
separation is achieved very rapidly conpn,red to the rate at 
which the electron density changes. As a result charge 
separation can always assxined to be in a steady state with 
no relative velocity between electrons and ions. This leads to 
the approximation that on a ma.croscopic scale the plasma is 
electrically neutral at every point and the electron density 
changes exactly equal ion density changes. In the gas dynamic 
shock the electrons, if present, also get compressed because 
of the strong electric forces exerted by the ions. With no 
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relatii/e -velocity, the electron pressure gradient arising 
from this compression cannot he supported hy collisions but must 
be balanced by the electric field. Thus, the body force on the 
electron due to the electric field must equal the electron 
pressure gradient. 

Together v/ith the regions of increased electron 
concentration in the front pa^rt of the shock front, there 
exists a region of reduced electron concentration in the 
rear part. This is because, the concentration of electrons 
without a simultaneous concentration of positive ions in a 
neighbouring region would lead to the appearance of an 
electric field at ’infinity* requiring the expenditure of 
an infinite amount of energy. 

1.9 St udy of Ionizing Shock Waves : 

r^Q Cn rjrj 

Earlier invest igators"^^^ ’ ’ ' have studied the 

relaxa.tion zone taking the initial conditions for the atom 
gas obtained from Rankine-Hugoniot relations for the inert 
gas shock. However, the choice of initial values for the 
electron and ion gases are not known a priori. Eor low Mach 
numbers when little ionization occurs inside the gas dynamic 
shock, the downstream equilirbium state of it is determined 
by the Rankine-Hugoniot relations and the initial electron 
and ion densities are zero. However, an initial electron 
temperature is still required. 



19 


The first in'vest igation of this type was carried out 
62 

by Petschek and Byron who ignored the effects of viscosity, 
heat conduction and radiation in their analysis. Collisional 
ionization was assumed to occur as a two-step process, the atom 
being excited and then ionized in a second collision. They 
considered only atom-electron collisional ionization and 
neglected recombination. The source of priming electrons was 
not clear from the work of Petschek and Byron. They used a 
steady state electron energy equation in which the energy lost 
by electrons in inelastic collisions is exactly balanced by 
the energy gained in elastic collisions with heavy particles. 

The energy transfer due to temperature gradients was neglected. 

Their calculation of the ionization rate was verified 
approximately by two independent experimental measurements. The 
first was a measurement of the electrostatic potentials due to 
electron diffusion in the gas. The maximum rate of change of 
the diffusion potential was measured with a probe and was found 
to be in reasonable agreement with the calculated rate. 

Continuum intensity measurements also served to check the 
calculated ionization rate. Radiation frora the recombination 
of ions and electrons was measured using a spectrometer and 
photocell with a fast response time. They also measured the 
time to reach equilibrium ionization. 

14 

Bond has also considered theoretically the rate of 
approach to equilibrium ionization behind shock waves in argon 
and calculated the ionization time for slightly stronger shocks. 
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He assumed "bhat atom-adom colxisions provide tbe initial supply 
of electrons. His estimate of the electron-atom ionization rate 
is greater thaa that of Petschek and Byron hut their results for 
the total ionization time are in agreement. 

Pollowing the work of Petschek and Byron, the 

experiments by Harwell and Jahn^^ and also Kelly^^ confirmed 

that at the Mach numbers considered the initial ionization 

mechanism was atom-atom collisional ionization. Morgan and 
50 

Morrison essentially repeated the work of Petschek and %-rou 
but included the atom-atom collisional ionization. 

¥ong and Bershader'^'^ made a similar analysis and 
included atom-atom collisional ionisation. They attempted to 
identify the dependence of the cross section on energy, for 
this mechanism using the available experimentally measured 
ionization rates. They also considered the steady state 
electron energy equation of the Petschek and Byron’s type. 

They used a slightly different expression for the atom-electron 
collisional excitation cross section from that of Petschek 
and Byron and included the three body ion-electron-electron 
collisional recombination. Their theoretical analysis was 
supplemented by experiments using a Mach-Zehnder interferometer 
in argon. Within the realm of their experiments 
photoionization, precursor effects and impurities did not 
seem to play a significant role. Recombination was found 
to play a minor role in deteimining the electron density 
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profile. Por Mach numbers greater than 18, the measured 
electron concentrations did not reach the Saha equilibrium 
concentration and radiation was stispected to play a significant 
role in cooling the gas in this region. 

36 

Ho.ffeif and Lien studied the relaxation zone by 
considering both atom-atom and atom-electron collisional 
ionization. The corresponding three body collisional 
recombinations were also included. They hawe used a complete 
electron energy equation which included the effect of electron 
temperature gradient. Por an initial vaJ-ue for the electron 
temperature they chose two values, one equal to the upstream 
oqLiilibrium atom temperature and the second equal to the atom 
temperature downstream of the atom-atom shock. Their nmerical 
solutions showed that only the early part of the relaxation zone 
was affected by the choice of initial electron temperature, and 
that the electron temperature varies veiy little in the later 
regions so that a steady state electron energy equation was 
sufficient. 

1 9 

Chubb studied the problem for larger Mach numbers 
(M'-^30) when ionization within the gas dynamic shock becomes 
important. In other words the width of the relaxation zone 
becomes comparable to the width of the gas dynamic shock. 

His approach is a semi-kinetic theoretical and the regions of 
interest namely the gas dynamic shock and the relaxation zone 
are studied separately. The gas dynamic shock is described 
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by a Mott-Smitb’ s type of solution for atoms and included 
a two-step atom-atom and atom-electron collisional ionization 
in both the regions. The corresponding collisional 
recombinations have been included only in the relaxation zone. 
He used separate electron continuity 5 moment-um and energy 
equations in addition to the plasma continuity equation and 
a c^ - moment equation^"^ for atoms. The plasma conservation 
equations of momentam and energy are not satisfied throughout 
the atom-atom shock| however, they were used to obtain jump 
conditions across the atom-atom shock thus satisfying 
momentum and energy conservation for plasma only across the 
atom-atom shock. In this case the quantities appearing in 
the structure equations for the atom-atom shock depend upon 
the electron density and temperature downstream of the atom- 
atom shock. Since these a2?e not known a priori for starting 
the numerical integration of the structure equations, an 
iterative method wa ,8 used. In the study of the relaxation 
Zone all the three species are a^ssuined uo have Maxwellian 
distribution functions. The three conservation equations 
for the plasma and a more complete electron energy equation 
than that of Hoffert and lien were used in which the velocity 
gradient tem was included. The set of equations were solved 
ntnnerically using for the initial values the downstream 
quantities obtained from the atom-atom shock. 



23 


As far as the gas dynamic shock is concerned for 
Mach nutmbers >20, 

i) the continuum equations cannot meaningfully 
describe the flow and a kinetic theoretical 
description has to be used; 

ii) the degree of ionization is large and inelastic 
collisional effects are to be included in addition 
to elastic collisional dissipative mechanisms. 

Por the relaxation zone earlier investigators have 
assumed that 

iii) each species may be described by its own Maxwellian 
distribution. This appears reasonable as long as 
the average thermal energies are low compared to 
excitation energies, and the density is relatively 
high so that, a particle will experience a large 
number of elastic collisions between the inelastic 
ones which change its internal state, Por the 
temperatures encountered in the relaxation zone 
behind shock v/aves for M > 30, and at low 
pressures the Maxwell distribution postulate 
for atoms is doubtful. 

1,10 Scope of Present Study : 

Chubb considered the aspects (i) and (ii) of 
section (1.9) hut did not satisfy the plasma momentum and 
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energy conservation throughout the atom-atom shock. In the 

present investigation not only all the above three aspects are 

considered hut also the plasma energy conservation explicitly 

satisfied throughout the entire shock finnt, further the 

iterative proceedure of Ohuhh is unnecessary in the present work. 

In addition assuming a nonequilihrium distribution function for 

atoms takes into account the elastic collisional dissipative 

mechanisms active towards the end of the relaxation zone 

where gradients in velocity and temperature are not negligible. 

The atom-atom inelastic collisional dissipative terms that 

have been neglected in comparison to the elastic collisional 

2 

dissipative terms in the c - moment equation of Chubb are 
likely to introduce considerable error at the end of the 
atom-atom shock. These terms have been included in the present 
study, A corrected value for the atom-atom collisional 
excitation cross section is also used, 

furthermore, entire region of nonequilibrium is 
represented by n single set of equations. This study may also 
be looked upon as an illustration of the suitability of a 
trimodal distribution function for a problem characterized 

by three equilibrium zones (the three zones do riot really exist 
in a physical problem). The radiation effects, precursor 
effects and multiple ionization have not been considered which 
are, however, important at higher Mach numbers and lower initial 
pressures. 
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COILISIOR PROCESSES AMD COLLISIOH CROSS SECTIONS 

2.1 Introducti on; 

Ionization is a process in which sufficient energy 
is imparted to the electron in its orbit around the nucleus 
so as to leave the orbit and pass into the continuum spectrum. 
This energy is known as the Ionization potential of the atom. 
Ionization may take place in mmre than one step in which case 
the atom passes through various excited states. If sufficient 
energy is available, each of the eleraentamy processes resulting 
in excitation of electrons in atoms can also result in 
ionization. Ionization is essentially a limiting case of 
electronic excitation. 

The elementary processes which result in excitation 
or ionization of atoms can be broadly classified into two 
categories; 

i) collisional process in which the required 
energy is transferred during the collisioH(s) 
between the atom and a high energy particle, 

ii) photoprocess in which the role of the high 
energy particle is played by a photon. 
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SEyniboli^ally the two processes can he written as 

follows, 

A + B a'^+ B + e 

+ ( 2 . 1 ) 

A + hv 5 =^ A + e 

where A is an atom , e an electron and A an ion. 

B is the secondary particle which participates in the collision. 
This may he an electron, which is so most often, another 
atom, an ion, or a foreign impurity particle (if the gas 
is impure) and hv is a photon. The in-verse processes 
represented hy the re-versed arrows describe the corresponding 
recombination processes. 

Strictly speaking, all these processes occur 
simultaneously in a gas. However, in certain situations, some 
of these processes may turn out to he less significant. In 
general, the process may follow a more complex path than the 
one-step process mentioned in eq.. (2.1), in which case an 
intermediate excited state of the atom also plays an important 
role. Such process can he symbolically represented as follows; 

A t A A*+ e + A 

* 

A 4- e A + e 

A 4* hv A* 


( 2 . 2 ) 
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and 

A + A A + e + A 

* + 

A + e A + e + e (2.3) 

* + 

A + tiv ~ — ^ A + e 

where A is one of the excited states of the atom, The 
reactions represented hy eq, (2,2) are those in vdiich the 
atom is excited to one of its possible lewels from the gro-und 
state, whereas those described hy eq, (2.3) represent the 
ionization of an excited atom, itirther, these processes may 
take place in multisteps. However, it will he shown later 
that the two-step process is the dominant one ever the rest 
at such temperatures that kT^I. In what follows, the 
atom-atom and atom-electron two-step collisional processes 
are considered in detail and the relative importance of these 
processes over the single-step process and also the 
photoprocesses is established. 

2. 2 Atom-Atom Oollisional Ionization : 

This is one of the dominant mechanisms during the 
initial stages of ionization when the electron concentration 
will be very low. In the absence of other initiating 
mechanisms like photo ionization, this mechanism determines 
the rate of increase of electron density in the gas^^. 
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Johnston and Komegay"'^ , employing microwave 

techniques demonstrated for xenon that atom-atom collisional 

ionization proceeds hy a two-step process in which the rate is 

governed hy the first process of the atom getting excited to 

its first level from the ground state. The rate of excitation 

c r ionization of the atoms depends on the number density of the 

particles that participate in the inelastic collision, the 

concerned collision cross section and the activation energy 

of the process. Harwell and Jahn^^ and recently McLaren and 
48 

Hohson have shown that their theoretically predicted results 
of the two-step model agree qualitatively with the measured 
ionization rates in their experiments. Assuming Maxwellian 
distribution for atom velocities, it can be shown that the rate 
at which atoms are ionized in unit volume through a one-step 
process is 



where n^, m^ and are respectively atomic particle density, mass 
and temperature j G the center-of-mass speed of the colliding pair 
of atoms and g their relative speed; g = ir(4I/m ); (jnCg) is 
the ionization cross section. Subscript 1 stands for one- step 
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3^0 0 6 3 3 • *4Jl ssuining a lineax* dependence of -bhe cross section on 
energy in the center-of-nass system (Section 2.5 ) as 

CTj = {E-I)H I (g^-g^), (2.5) 

eq.. (2.4) can he vc^itten, a.fter carrying ent the integration;, as 



1 



whe re C 


(a) 


2 

B. 3. 


is a constant and a is the rate constant. 
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of the temperature. 


I 



( 2 . 6 ) 


(2. 6a) 
a function 


The rate of excitation of atoras, per unit -volune from 
their ground state to the first possible energy le-vel can be 
written in a similar way, by replacing in eq. (2,6) 0^®'^ and I 

( Sl) 

by the corresponding constants C4 ' and as 
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1 + 


a^ 
(2.7) 




In expression for the rate o.f ionization )f the excited atoms 
can be similarly written o.s 
dn , 


■ST 


= 80 


(a) 


in n )■ 

C.I 


■V2 


1 + 


2kT 

a 

( 2 . 8 ) 


I’ 


where 0^^^ is the constant of proportionality in the cross 

T* 

section for ionization of the excited atoms, (E-I') 

and I* = I-I* is the activation energy for the process. 
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Since for the noble gases even the lowest excited states 
have energies well above half the ionization potentials, the 
energy increments involved in the second step, I' ’vill be subs- 
tantially less than those for the first step, 1.,^., Also the 
geometric size of an excited noble gas atom presumably exceeds 
that of its ground state, and a continuura of final states is 
available for the second step compared to a discrete level for 
the first, for this reason we can assujme the chances of an 
excited atom getting ionized in a further collision to be large 
compared to the probability of collisional deexcitation. The 
process of ionization of an excited atom is faster than the 
process of excitation of an atom from the ground state. 

The rate of ionization of excited atoms hinges on the 
steady state concentration of the excited atoms A (eq. 2.8). 
This, in turn depends on the relative rates of excitation and 

deexcitation processes. The natural life-time of a noimal 

— R 

excited atom, against radiative decay (>^10 sec) is perhaps 
inadequate to support the required steady state concentration 
of A*. However, in noble gases, one of the following aspects 
may effectively protract this life-time. In noble gases like 
argon at moderate pressures and temperatures the excitation 
collision results in the elevation of the atom from* the ground 
state to one of the four closely spaced energy levels comprising 
the first excited state. Two of these are metastable states which 
are energitically isolated from the ground state and any 
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transition between these levels and ground state are forbidden 

by virtue of the J selection rules. Spectroscopic transitions 

between the other two levels and the ground state are permitted 

by the selection rules and give rise to so called resonance 

radiation. If collisional excitation leads to one of these 

metastable states, the atom can remain in that state for periods 

of the order 10”^ sec.^^, because of its inability to reach the 

ground state by radiative decay. Consequently the chances of a 

second collision leading to ionization of this atom is greatly 

enhanced. If the atoms are collisioiially excited to one of the 

resonance levels, which, for the noble gases studied lie less 

than 0.3 eY above the lowest metastablo states, the life-time 

of the excited atoms for radiative decay is of the order of 
8 

10“ sec., which is usually less than the mean collisional time. 
However, since the gas is optically thick for this resonance 
radiation, the emitted photon will be trapped in the gas, indi- 
vidual atoms being alternately excited collisionally or radia- 
tively and deexcited radiative ly resulting in a constant over- 
all population of the excited atoms. Therefore, the rate of 
ionization of an excited atom by a collision of the second type, 
which depends on the population of n*, increases with resonance 
radiation trapping. The role played by this resonance radiation 
is discussed in Section (2.42). Since the difference in the 

4 

energy levels between the metastable and resonance levels is 
commensurate with the experimental inaccuracies it has not been 
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possible to experimentally verify which of the tvro paths the 
mechanism really follov\rs« However, the experimental results of 
McLaren and Hobson^ suggest that met 8.st able states play the 
main role in the mechanism. This however needs further inves- 
tigation. Because the activation energy for the forv/ard reaction 
leading to ionization of excited atoms is considerably less than 
the activation energy for the deexcitation process, the ionization 
is more probable than the deexcitation of excited atoms. Due to 
the large difference between the activation energies of the exci- 
tation reaction and the ionization reaction of tho second type, 
the former being several tines the latter, the former process 
controls the net rate of production of electrons. 

Since the energy req,uired t o excite an atom from the ground 
levels is less than that required to ionize an atom in one- step 
(l^ < I) and because the number of particles that have energies 
above is much large than those ha.ving energies greater than I, 
the number of collisions that result in tho excitation of atoms 
is large compared to the number that result in direct one-step 
ionization. Further, almost all the excited atoms get ionized 
by a second collision. Thus the two-step ionization process is 
the dominant mechanism in noble gases and the excitation colli- 
sion process governs the overall rate of ionization. 

Experimental measurements of the ionization rates in 

53 3Q 30 

noble gases by Harwell and Jahn-^-^, Kelly-'^-^, Margon and Morrison 

/ Q 

and WicLaren and Hobson^ also suggest the choice of the two -step 
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ionization mecbanisni. Following fhe analysis of Harwell and 
33 

Jahn 5 for a noble gas free from impurities, the possible 
initial two-step ionization mechanisms can be described sym- 
bolically as 


A + A 


A + A 


A + hv 


(2.9) 


+ A 


e + A 


4. 

where k*, k*' , k^ and k are the corresponding reaction rates, 
which are functions of the temperature; (2,9) represents the 
photoexcitation and deexcitation mechanisms. 

Taking natural logarithm, and differentiating the eq,. (2.6) 
with respect to (1/kT) yields, for kT << I, 


d(1/kT) 


( 2 . 10 ) 


expressing the usual Arrhenius result for one-step process. 


If the ionization takes place follov/ing a two-step path 


then. 


dn , 

( tA ) 


dt ^ 2 dt 


( 2 . 11 ) 


Again, from eq. (2,7) we can write the expression for the slope 
of the Arrhenius plot as, 


d(l/kT) 


( 2 . 12 ) 
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Measurements have shown the slopes of the In "^s. 

(l/kT) curves to he closer to than I, which confirms the 
relative importance of the activation energy over I. lor 
example, in the measurements made by Harwell and Jahn using 
a transverse microwave detection arrangement with Ar, Kr and 
Xe, the activation energies deduced from the Arrhenius plots 
are 11,9 + 0.5 eY, 10,4 + 0.5 eY, and 8.6 + 0.5 ©Y respectively, 
which are closer to the corresponding excitation energies, namely, 
11.55 to 11.83 eY, 9.91 to 10.6 eY , and 8.31 to 9.55 eY for the 
respective gases; these activation energies differ distinctly 
from the respective ionization potentials 15.75 e\, 14.00 ey, 
and 12,13 eY. Recent experimental investigations by McLaren 
and Hobson^ using double electrostatic probes in argon have 
revealed further aspects of the two-step process, fhe activation 
energy got from the Arrhenius plot has a value 11.4 + 0.3 eY, 
which is closer to the metastable level energy of o.rgon 11.6 eY. 
They were also able to measure the activation energy I' for ioni- 
zation of the excited atoms and the corresponding rate constant 
This activation energy was found to have a value of 
4.7 + 0.4 eY, which compares closely with 4.27 eY, required to 
ionize an excited argon atom from its metastable Mevel. These 
observations augment the model of a two-step ionization mechanism 
involving an intemediate metastable state of excitation. Using 
the values of the mieasured reaction rates k’'"’ and k"^ they were 
able to estimate the steady state degree of excitation 
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n^/na (= k^/k^) over the range of tempera,-tures from TOOO^IC 

o ^ 

to 12000 K and at pressures from 0.2 to 5.0 torr. The concen- 
tration of excited atoms (n.-Vn ) was found to vary from 7 x 10~^ 

^ 3/ 

to 8,64 x; 10“ vmicb is sufficient to keep the rate of ionization 
of excited atoms much greater than the ra,te of excitation. Ratio 
of the cross sections for the two processes found 

3 

to be v,'' 10 for argon which is of the same order as estimated by 

58 

Johnston and Komegay for xenon . 

55 

Harwell and Jahn"^ also highlight the relative importance 

of various deexcitation processes by studying the dependence of 

reaction rates on time and density. The rate equations for the 

excitation and ionization processes represented by (2,9) can be 
55 

written as 

2 - 
U- = k* na - (k^ + k^’) n^n^ - k^ n* 

(2.13) 

dn^ 

Assuming n^ to be constant which is a good approximation 
for veiy lovif degrees of ionization, equations (2.13) can be 
integrated assuming n^ , n_^ = 0 at t = 0 to give, 

= k^ k^T na^ [ t-T(1-e“‘^/'^)J 


where x 


(k^ + k*') ^a + k^ 


(2.14a) 

(2.14b) 

(2.14c) 
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The results reproduced here for ready reference fron (33) 
and shown in Table 2,1. give the activation energy, the density 
and time dependence of the electron concentration deduced from 
equations (2.14) for different possible cases. 

Their experiments havo shovm that the density dependence 
waa quadratic and the time dependence viras linear, though the 
uncertainty about quadratic dependence on density was expressed 
to a certain extent. In some cases they found density dependence 
weaker than quadratic, which any way favours none of the alter- 
native paths other than the second shown in the Table (2.1). So 
within the realm of the investigations of Harwell and Jahn, it is 
justifiable to neglect the collisional and radiative deexcitation 
processes. 

i'iCLaren and Hobson' a.lso have studied the relative import- 
ance of collisional deexcitation over collisional forward ioniza- 
tion of excited atoms. They show? from the consideration of 
energies involved in the tv/o processes, that neglect of colli- 
sional deexcitation compared to forward ionization requires that. 


>> 




(l^kT) 


(2.15) 


“0 ■"+ 

where g^ and g^ are the degeneracies of the ground and excited 
states. In the range of their experiments in argon it is shown 
that the quantity on the right hand side of the inequality (2.15) 
varies from 0.05 at 12000°K to 1.2 at 7000°K. For the metastable 
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TalDle 2.1 


Various Ionization Routes in the Initial Regime 


Process CollisionaR Rad. decay Phase Density Time Activa- 

De-excita~ Order Dep. tion 

tion energy 
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state of the argon atou the ratio (g^/g^) has a minimm value 
of 6 which again justifies the neglect of collisional deexcita- 
tion process. 

2 . 3 Atom-Electron Collisio na l Ionization : 

In this process high energy electrons collide with neutral 
atoms res-ulting in transfer of energy inelast ically to the atom 
sufficient enough to excite the atom to its first excited state. 
Because of the fact that the mass of ai" electron is several 
orders of magnitude less thaji the mass of an atom and that the 
energy transfer by collisions depends upon the mass ratio of the 
collision participants, the efficiency of this process is very 
much different from the atom-atom collisional process. 

i'he condition for the probability of inelastic energy 
transfer during a collision to be high is that the adiabatic 
factor (27 t ; au/v) should not be too great' conipared to 1 , where 
•a’ is Bohr radius, m the angailar frequency of an electron in 
orbit, and v the relative velocity of the approaching particles. 
Thus, in order that a collision be inelastic the impact must be 
sufficiently sharp, i.e. the velocities of the approaching par- 
ticles must be of the same order as the velocity of the electrons 
in the outer orbit of the atom. In the case of an inelastic 
collision between two atoms -with relative energy of the order 
of the ionization potential, say upto 10 e7, this condition is 
not satisfied and the inelastic cross section is very small. On 
the other hand when an electron approonhes an atom with relative 
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energy of the order of I, tiae -velocity of the ele ctron5'\r(2I/m^ ) , 
will he of comparahle order to the -velocity of the electix^n in 
orhit thereby uaking this process efficient in ionization. 


Expressions similar to (2,6) can be written for the rate 

of generation of electrons by the aton-elc-ctron collisional 

process assuming a liner-ar variation of the cross section with 

62 

respect to the energy in the relative ^stera as 


dn 


(e) 


am -1/2 


4 Cj ^ ) 






) e 


• I/kl' 


^‘e ^a^e 


fe 

(2.16) 

(2.16a) 


.(e) , 


where Cj is the proportionality constant for the cross section 
given by, 

-(i) ^ Q(e) 


(T. 


ae 


(E - I) 


(2, 17) 


where n^, are the mass, density and temperature of the 
electrons; is the rate constant, a function of the electron 
temperature. 


Comparing a_ and from (2.16) and (2.6) for the same 

V- 3. 

electron and atom temperatures. 


a. 


a. 




c(e) 


/ m n^ 
/ _a _e 

,/ m^ n 
e a 


■f(m /m. ) 100; and measured values 

^ 10 ^. 


(2.18) 


of C 


('^) 


and 


p ( ^) 


give 
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Hence as long as n^/n^ is not too small (not 

less than 10~^ - 10“^). Liorgan ct al,^° have foimd that the 

take over degree of ionization varies from about 5 x 10“"^ for 

temperatures of 10000°K to 5 :r 10“^ at 40000°K;, Vdien 1 is 

e 

much less than atom temperaturesj the onset degree of ionization 
takes values greater than these. 

I'or ionization fractions greater than 10"^, the electron- 
atom collisional ionization process dominates over the atom- 
atom process. In most casesj the initial building up of this 
electron concentration takes place through one of the priming 
processes. After the take over point, the electron-atom colli- 
sional process becomes dominant over the rest, 

Ihe electron-atom collisional ionization also takes 
place in two steps similar to atom-atom process, involving an 
intermediate excitation stage and the rate of ionization is 
again governed by the rate of excitation of the atoms from 
the ground state. The ratio of the characteristic time for 
ionization to that of excitation of atoms by electron impact 
for argon can be written as'^® 


xe 

T* 

e 


1 

10 



el/kT 


for l/kl V/-' 10. 


( 2 . 19 ) 


Close to equilibrium and with moderate degrees of ioni- 


zation. 
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11 . 


e . 


n 


( 


a 


a 

For argon at n = 1.7 x 10 
a 


^ 3/4 

> -^eT ' 


n 


( 2 . 20 ) 


S *5 
CEi"''^ and 


135000°K, 


'^ q/'^q '■^ 5000 indicating that the Boltzmann distribution with 
respect to excitation is established much faster than in ioniza- 
tion, The times may be comparable only at the beginning of 
the process when the electron density is much smaller than the 
equilibrium -value. This along with an argument similar to that 
followed in Section (2.2) shows that the excitation rate is 
much larger than the ionization rate and that the excitation 
rate dictates the ionization rate in the two-step process. 


2 ^ Pho'bo ioniza t ion : 

Photoprocesses become of considerable importance at high 
temperatures in determining the shock structure. The radiation 
emitted at the end of the relaxation zone as a result of ion- 
electron recombination j propagates upstream towards the cold 
gas. In the process it pairtly or completely gets absorbed by 
the gas. This escape of radiafion from the downstream equili- 
bri-um region and its partial absorption by the gas results in 
(i) cooling of the downstream equilibrium gas and (ii) change 
of the shockwave structure due to photoexcitation and ionization 
of atoms. 


The photoprocesses that result in the emission of radia- 
t ion are 

(i) free-free transitions or Bremsstrahlung radiation. 
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arising iron ihe deceleration of the electron in free state 
in the conlomh field of the ion with the emission of a photon 

+ e + e + hv (2.21a) 

(ii) free-ho-und transitions, which inwolre the photo recombina- 
tion of the free electrons into one of the bound states of the 
atom, the difference in energy being releasedsc as a photon 

a"*" + e “ a(p) + hv (2.21b) 

where A(p) is an ezcited atom in p-th electronic state, and 

(iii) boimd-bo-und atomic transitions, in which an electron 
from a higher energy level moves to one of lower energy wfith 
the emission of a photon 

A(p) - A(q) + hv (2.21c) 

Since the electrons in the free state can have a continu- 
ous distribution of energj^ the first two processes v^hich involve 
transitions from free electronic state, result in continuum 
radiation while the bound-bound transition involves discrete 
energies and hence emits line radiation. As far as the energy?- 
loss and cooling effect are concerned, the continuum radiation 
maAes a large contribution. Various aspects of the emission 
and propagation processes and the resulting cooling effects of 
plasma have been considered by several authors in References 20, 
27,34,37,47557,59,63,64,67,70,72. In this section, the relative 
importance of photo ionization vis-a-vis collisional ionization, 
is considered and the resulting cooling effect has been ignored. 
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The former effect influences the entire shock structure whereas 
the latter effect only the tail cud regions, 

Effect of Continuum Radiation; 

Assuming a one- step photoionizaf ion , the rate of produc- 
tion of electrons due to photo ionization can be written as'^'^, 

CO . . 

Z_ = n r ^^-4 c cr(v) dv (2.22) 

a. J hv V ^ ' 

where U(v) is the radiant energy density per unit spectral 
interval at v , (J^ (v) is the phato ionization cross-section for 
the atom from the ground state, and c is the velocity of light. 
Since the energy of the photon must be greater than I to ionize 
an atom, the lower limit for v is given by hv^ = i. 

_3 

has a frequency dependence of v for hydrogen-like atoms. In 
complex many-electron atoms the assumption of hydrogen-like 
behaviour for excited states of large principal quantum ntanbers 
is reasonable, but the ground state and low lying excited states 

O 

are highly nonhydrogenic. Bates and Seaton , through detailed 
quantum mechanical calculations, have established that the 

3 

frequency dependence of can differ markedly from v depen- 
dency for ground state photoelectric effect because of the 
non coulombic nature of the field between the photon absorbing 

electron and charges due to nucleus and the remaining electrons. 

9-11 

Bibor'ian ot al. developed approximate analytical expressions 

for free-free and bound-free transitions from excited states, 
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for a variety of low atomic number atoms* I'be expressions 
are valid for local thermo dynamic equilibrium and resemble 
the well known Kramer-Unsold expression for hydrogen. The 
effective nuclear charge Z has been repla,ced by a correction 
factor ^(v s T), that depends in general, upon both fre- 
quency and temperature. The factor ? (v , T) has been evaluated 

1 6 

by Bibeiman and Norman using Burgess- Seat on theorj’’ , for a 
number of neutral and ionized atoms. With the crudeness of the 
following analysis, which is intended only to estiiiiate the 
order of magnitude of the photo ionization rate, an average 
cross section (which can be set equal to the cross section at 
the ionization threshold to a high degree of accuracy) can be 
used. The Planck equilibrium radiation energy density can be 
used for U( 

U(v) = (e^^' (2.23) 

c^ 

further, since (l/kT)>> 1, (2.23) can be approximated by the 
YiTien ’ s di splacement exp re s sion , 

0(v) g-hv/lcT 

C^ 

Using Cr^(^) = a^nd substituting from eq. (2.24) for 

U(v) eq. (2.22) can be written for as, 

Bren KT-y 2 P 2 ^ T /kT 

Z . (Tv. — ^ (I + 2IkT^ + Bk^T*^ ) e 


(2.25) 
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where is the eq.uilihriim temperature for the radiation. 

Putting I >> kP simplifies eq. (2.25) to. 


l2 -I/kT„ 


(2.26) 


in expression similar to eq, (2.16) for the two-step atom- 
electron collisional ionization rate can he written for I^>>kTg, 


/ ^ s 8kT^ 1 /2 

.(e) / e X ' 


■Ik At, 


A- ^ I* e ' ' ® (2.27) 


a e 


where ' is the constant of the cross section expression for 
the atom-electron collisional excitation. Prom equations (2.26) 
and (2.27), for argon^ using the values = 11,5 eV and I = 15.7 eV, 


^ ... 1-18 X exp 

h “e iRtg ^(e) 


(2.28) 

he - 


In the regions where electron-atom collisional ionization domi- 
nates over other collisional processes, Tg T^^^ (= T^). There- 
fore eq. (2.28) can be written as. 


'^^1 1.18 X lo"^^ ^V1 


XeT ^'Aq 


fAn ®^P ((- ^”)X10 

eq 


(2.28a) 


Por an equilibrium temperature of 10^ °K (for an upstream Mach 
number of 30, and temperature 300°K, the equilibrium temperature 
is 17100°K) the ratio of photo ionization to collisional ioniza- 
tion is. 


^vi 2 X lo"*^ ^1 


3F ■ 


(2.29) 
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Ass-uming ^ fTovn eq. (2.29) it 

is foimd, that photoionization is relatively less important 
compared to collisional processes. For lower electron densi- 
ties however, the photoprocesses need to he considered. 

At high Mach numbers (LI >20) there exists still quite 

an amount of -uncertainty about the significant mechanism of 

ionization in the initial stages which build up the electron 

population to the onset value. For lov/er Mach numbers, it has 

been fairly well established^^ that atom-atom collisional 

1 Q 

process dominates over the rest. Chubb estimated the ratio 
of photo ionization rate to collisional ionization rate assuming 
a constant photoionization cross-section and equilibrium radia- 
tion corresponding to the downstream equilibrium. His analysis 
is similar to the one followed above and for ^ he has used 
the largest value of photoexcitation cross section quoted by 
lee and Weissler^ for argon. Chubb's expressions and calcula- 
tions, which seem to be in arithraeticsA error* have under esti- 
mated the significance of photo ionization. The corrected value 
of the ratio of photo ionization rate to atom-atom collisional 
ionization rate for a Mach number of 20 is found to be 5.2, 
thus it seems that photoionization is important even for Mach 
numbers around 20 for the dciSgit-fus considu-r^I-d, However 


* Ref. 26, Eq. 50 is less by a factor 10 1 and the value of 
Ip/lcoii M = 30, quoted as 5 x 1 0~2 is veiy much off the 
correct value of 11.2, 
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such a conclusion cannot be made on the basis of the crude ana- 
lysis of the above type. The assumptions of constant photo- 
ionization cross section and eq.uilibri-u]ii distribution of radiant 
energy are highly questionable, and the lowering of equilibriira 
temperature because of radiation cooling has not been considered. 
Since the cross section decreases for frequencies above the 
threshold value the constancy assimption overestimates the 
photoionization rate. Secondly, the effect of radiation absorp- 
tion by the nonequilibrium gas as it passes through, before the 
radiation reaching the atom-atom zone is not considered. Kelson 
and Goulard"^' have considered the effect of absorption of the 
ground state continuum radiation for argon gas for initial 
pressures of 0.2, 1, and 10 cins of Hg and an initial temperature 
of 300°1C for a Mach number of 24. In their analysis the 
electron gas was allowed to be colder with respect to atom gas. 
Within the realm of their studies it is found that absorption 
effects are important to the extent that ground state continuum 
radiation does not affect the region behind the shock much and 
the zone is collision dominated. The absorption effect can be 
studied by considering the non dimensional optical thickness 

o 

nondimensionalised v\rith respect to collisional ionization mean 
free path, which they have calculated through-out the shock. 

When r„ is large, the ra,diation emitted by a particle is 

g 

essentially absorbed by its neighbouring particle, so that 
the particles just exchange equal amounts of energy and no 
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net loss or gain occurs; consequently, the radiation is trapped. 
As the optical thickness decreases belov; unity, the radiative 
energy can travel quite far before being absorbed. The results 
of ITelson and Goulard shov/ that x > 1 in the atom-electron 

o 

dominating region and Tq, < 1 in the atom-atom dominating 
region, so that the groumd state continuum radiation emitted 
in the atom-electron region is practically trapped and only 
that emitted very close to atom-atom region reaches the upstream 
region and participates in ionization. The radiative fluxes 
and optical thicknesses depend strongly on the electron tempera- 
ture; thus the cool electron gas in the atom-atom region greatly 
curtails the magnitude of the radiation entering the atom-atom 
region. The assumption of equilibrium distribution of radiant 
energy is invalid; because the relative changes in temperature 
and density a.Te large within distances equal to the penetration 
lengths pertinent to the radiation, nonequilibrium radiation 
may occur. It may also occur in optically thin regions, because 
the radiation emission rates are not balanced by the absorption 
rat<3 s ♦ 

2.4.2 Effect of Line Radiation ; 

The one-step radiative processes discussed above exclude 
the effect of resonance line radiation resulting from emission 
processes of the type (2,21c). This process has often been 
dismissed due to the strong self-absorption effect expected at 
the center of the line. Nevertheless, the lov/ absorption 
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characteristics of the wings of the line give rise to the 
possibility of a substantial optical path for the escaping 
photons of the sane frequency j until they are trapped? thus 
in the atom-atom region and also in the precursor region, the 
excited states formed by the absorption of this line radiation 
may get further ionized by the absorption of a second photon 
of continuum excited state radiation. 

12 *51 

Biberman and Veklenko and Murty- have theoretically 
investigated the importance of this two-step precursor mechanism, 
d'he experiments of Wilson and lederman"^^ in argon showed stronger 
precursors at higher densities v/hich is the reverse of what any 
one-step ground state continuum theory would predict. Dobbins' 
two-step model^^ and also Murty 's^^s on "the other hand depict 
the same experimental trends, though q.uantttatively they are 
much far off. In their results the ground state continuum 
curves are also given which suggest that both one-step and 
two-step processes coexist atleast in argon gases. However, 
much work, particularly on the experimental side needs to be 
done to understand the radiation exchange process in argon 
and other gases. In the absence of experimentally measured 
photoionization cross sections, the atom-atom collisional 
ionization process is considered as the dominant initiating 
mechanism in the present investigation. 
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2 • 5 I nelastic Col lis ion Gr oss Sec t J.o_n s i 


In tills section the inelastic collision cross sections 
lor the atom-electron and atoin-atom collisional excitation 
will he considered. The first atteimpt to study the classical 
atom-electron inelastic collision cross section Vvas made hy 
Thomson in 1912. He described the inelastic collision on the 
basis of classical mechanics assuming thaf the collision of 
a free electron with an atom, takes pl3.ce in a time short 
compared to the period of rotation of a bound electron in its 
orbit so that the latter acquires energy from the impact in 
the same manner as would a free electron. If an electron with 
kinetic energy E passes near a target electron, the cross 
section for the energy transfer to the target in the range a E to 
A E + dA E from classical mechanics is. 


dO- = Esl (2.30) 

E (AE)^ 

where e is the electron charge. The cross section for the 
transfer of energy E, greater than a threshold value E^ is 
therefore , 


^4 .E 

cr= 2£|_ j 


A 


dAn Tce \ 

-p— I — — ^ • 

^ “ E (AE)^ E 


(2.31) 


Eor E not far from E^, eq. (2.31) can be written in 


the form, 


Cr^G (E-E^) 


( 2 . 32 ) 
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where C is a constant and eq_, (2.32) implies linear dependence 
of the reaction cross section on the energy in the center-of- 
mass coordinate. Most often the value of 0 has "been determined 

experimentally for various gases. Petschek and Byron^^ in their 

** 

analysis, determined the value of this for argon hy making 

a linear approximation to the inelastic cross section data of 

Maier and leihnitz^^, vjhere as V/ong and Bershader*^"^ used a two 

slope linear approximation. In the present work the value of 

fiP 

Oae reported hy Petschek and Byron has been used. 


Experimentally measured values of cross- sections for 
atom-atom collisional ionization are available in t he litera- 
ture"^ ’ ^ They assume a linear dependence of cross 
section on energy and attempt to determine the constant of 
proportionality using the measured reaction rates. Wong and 

rjrj 

Bershader attempted to identify the relationship between the 
cross section and energy. However, as pointed out later by 
McLaren and Hobson, they seem to have erroneously a.ssmed the 
constancy of a factor (Cg/p ) in their expression for the cross 
section (See eq.. (3) of Bef. ( 77) ) v/hcreas is, infact, 

foimd to be the constant term over their range of temperatures. 

With this modification incorporated, the cross section becomes 
temperature dependent and is, therefore, not of general applica- 
bility, Recently McLaren and Hobson^® have considered constant, 

linear and quadratic dependences of the cross section on the 

Sub scripts ae refer to atom-electron collisional prnc^s and 
aa refer to atom-atom collisional process. , 

.. ! . KANPUR, 

,jeNT«AjL uwhahYi 

! ‘54"r 


I 
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energy in the center-of-mass syatea. They have calculated the 

values of C for the three cross sections from the measured 
aa 

reaction rates for different temperatures. Their results show 
that for linear dependence, C varies least with t emperat-ure , 

If the cross section for which C changes least with tempera- 
ture is taken as the criterion, this establishes the credibi- 
lity of linear cross section. However, as the cross section 
vs. energy curve has a peak at very high energies after which 
the cross-section starts falling with increasing energy, the 
linear approximation will be in error at very high energies. 

In the present investigation a linear dependence of cross 
section with relative energy is used assuming that all the 
particles with energies just greater than I* result in an 
inelastic collision and lose their energy and thus will seldom 
acquire energies much greater than I^^. 

dd — 1 9 

Harwell and Jahn have reported a value of = 7x10 

p errs 

cm /eV for argon gas. Morgan and Morrison^^ carried out a theo- 
retical reassessment of the ionization mechanism and later 
referred to the experiments of Harwell and Jahn with some equip- 
ment modification, which showed that the reaction rates reported 
by Harwell and Jahn were too high. They showed that values of 

C should be reduced by a factor 10 to get the best fit rela- 

39 

xation times measured by Petschek and Byron. Later Kelly"^-^ 
from apparatus substantially the same as that of Harwell and 
Jahn, but with one order higher degree of purity of the gas 
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arrived at a value for = 1.2 x 10“^^ cn^/eV. McLaren and 
4-R 

Hobson^ in their experiments have incorporated corrections to 
allovv for flo?/ nonuniformities resulting from the effect of 
boundary layer growth on the shocked gas flow. As reported by 

McLaren and Hobson, had the appropriate corrections been made, 

20 

Morgan and Morrison’s result would come close to 2.5 x 10"" 
cm^/eV while that of Kelly to 4 x 10“^*^ both in good 

agreement with the value measured by McLaren and Hobson namely 
Caa = 2.5 x 10"^^ cm^eV. 

20 2 

In the present work the value C ^ = 2.5 x 10“ cm /eY 
has been used for argon. Since the values corrected for flow 
nonuniformities are not available elsewhere for other gases, 
the ratio of the corrected to uncorrected values of 0 „ is 

aci 

assimed to be the same for all the gases and the values reported 
by Kelly with the correction applied are used. Since the values 
of C „ for gases other than argon are not available elsewhere, 
it is assumed that the ratio C is the same as that of 

argon for all the gases to evaluate 0^^^ for other gases. This, 
however, may keep the quantitative results for gases other than 
Ar in error which of course will nut be severe due to the fact 
that the shock structure is veiy insensitive to the value of 
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2,6 Re c oia'b in at ion Ble chan i s n : 

7 

Bates, Kingston and licWhirter ' ba.ve studied the general 
recombination meclianisni in a plasma considering all the processes 
comprising electron capture, electron impact ionization, discrete 
transitions to distant le-vels, radiative capture and spontaneous 
radiative transitions. The collision processes predominant in 
a dense plasma and the radiative processes predominant in a 
tenuous plaaiia are two limiting cases of this general 'impact- 
radiative recombination’ process. The three body collisional 
recombination process studied by Bates et.al. ,can be used for 
gases like argon. In accordance with their model, the upper 
states tend to equilibrate rapidly with the unbound continuum 
states. Thus the net recombination rate equals the ground state 
atom formation by collisional deact iva.tion and spontaneous 
transitions from excited states as well as by three-body and 
radiative recombination processes from the continuimi. The prin- 
cipal difficulty in evaluating the imp act -radiative recombination 
rate by an appropriate model is that the cross sections and 
transition probabilities of these processes are unknown. However, 
in a dense plasma the recombination rate is controlled mainly 
by collisional transitions to the ground state. In the present 
analysis, only three-body recom.bination to ground state is con- 
sidered assuming all the excited states, in turn, to be in equi- 
librium, This recombination rate is related to the ionization 
rate through Saha equation. This assumes the concept of chemical 
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equilibriiHn 1:0 be preserved ir, the presence of -bbennal non- 

equilibriian. Strictly speaking this assimption is valid under 

chemical equilibrium but extension to nonequilibrium states 

are generally accepted. Further, as recombination is important 

only at the end of the relaxation zone v^/here the electron and 

plasma temperatures are not veiy different from the equilibrium 

temperature, the error arising from this assumption will not be 

large. Thus the ratio of the recombination rate constant to 

ionization rate constant for the atom-atom and atom-electron 

processes and can be written, 

0 / 0 . 0,0 

(2.33) 

= (2ZyZ^)(2TtmgkTg/h2)V2 gxp (-l/kT^) 

In the expressions (2.33), Z^ and Z^ are the electronic 
partition functions of the singly ionized ion and atom respec- 
tively which can be replaced by g and g„ the ground state 

"F V 

degeneracies of the ion and atom respectively without much 
error. Thus eqs. (2.33) reduce to, 

Kaa = (2g^/gg_)(27i:mgkT^/h2)V2 (-lA^^) 

(2.34) 

Kae = (2g^/g^)(2mmgkTgA^)^^^ exp (-I/kT^) 

2 • Elastic Collision --Oross Section s : 

2,7.1 Atom-Atom Elastic Interactions ; 


Several models have been proposed to represent the 
interaction between neutral atoms and molecules, A long-range 
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attractive and short-range repalsivc type of force la,ws with 
experimentally determined constants have been found to he 
quite appropriate in describing the ma-croscopic transport pro- 
perties of the gases. At high temperatures, what is required 
is precise infoimiation concerning the repulsive part of the 
interaction since high temperature properties are relatively 
insensitive to the attractive part, in inverse-power repulsive 
force law model would be adequate for the high temperatures 
encountered in shocloffaves. Expressions for interaction poten- 
tials obtained by fitting curves to the experimentally measured 
potentials at high temperatures are reported for several gases 
in references 2, 5j and 4.. Amdur and Mason-^ have computed 
the viscosity of various gases at high temperatures using the 

potential measurement curves of references 2, 3 and 4. Recent 

•1 

measurements of viscosity by Aeschliman and Cambel are fonnd 

to be 3 to 10 percent greater than the values reported by Amdur 

and Mason in the temperature range 3500° ~ 8500°K. Amdur and 

Mason also suggested that for high temperatures a more realistic 

model should include a factor in the potential. But they 

did not use it because collision integrals for the exponential 

potential had not been evaluated at that time. Later computations 
24 

of Levoto with exponential potentials using the integrals 
evaluated by Monchick"^^ show better agreement of viscosity with 
Aeschliman and Cambel' s values to within 5 percent. It thus 
appears that the exponential potential may indeed be a more 
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realistic representation of the collision process. However, 

the complexity of the integrals involved in the shock structure 

problem precludes the use of these realistic models so that 

often simple inverse power potential laws ha.ve to be resorted 

to. We use an inverse- fourth power law of potential in the 

present work choosing the value of the const ant in such a way 

that this law approximates the potential laws reported in 

references 2, 3, and 4 with least error in the range of interest. 

Thus the atom-a.tom force law is, 

B 

B = -^ (2.35) 

r^ 

where r is the distance between the two gas atoms during an 

encounter and B „ is the force law constant, 
aa 

2.7.2 Atom-ElectTOn El as tic Interactions : 


In order to describe the main features of the electron 
neutral atom interaction, the wave like np.uure of the electron 
should be introduced. This wave length called the de Broglie 
wave length is defined by, 



( 2 . 36 ) 


where v is the total classical velocity including any variation 


due to the field of the atom. When this v/avelength is comparable 


to atomic dimensions interference phenomena can occur that result 


in rapid variation of cross section with energy and angle. This 
effect is called the Ramsauer-Townsend effect, Eor low electron 
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temperatures, the long range polarization field can he used 
effectively. At electron energies for ?7hich Ramsauer-Tovmsend 
effect becomes important the interaction can no longer he re- 
presented hy a force lav/ of the polarization type, in which 
case there exists a Ramsauer minimum for the cross section. 

Several investigations both experimental and theoretical 
have been made, e.g. ref. 28,60,65 and 66 and the cross sections 
over a range of electron energies are available. However, these 
cross sections are not easily amenable for analytical handling 
in such complex problems as the present one, so that approxima- 
tions have to be resorted to. Since the electron-ion elastic 
energy transfer is dominant in heating the electron gas over 
the entire shock except, perhaps, in a narrow region in the 
beginning, (more so with the assumption of thermal eq.uilibrium 
between electron and atom gases ahead of the shock) the electron 
temperature profile is very insensitive to the choice of the 
atom-electron cross section. Hence, an inverse fourth power 
law of polarization field is assumed. It should also be noted 
that in the beginning of the shock where the atom-electron 
collisional energy transfer is dominant, the electron tempera- 
ture is low enough to justify neglecting interference effects. 
Thus, for argon at Te 300°K the momentum transfer cross section 
from the polarization law is = 8,1? which compares 

with the value 8,0 of Pack and Phelps^ (O'Malley ’ s^*^ results 
give a value of 14). At T^, v.- 1000°K, 

O'Malley's results give a value of *^1. 


= 1 . 425 vtoereas 
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Table 2.2 

Various Properties and Gross Secti ons for the G-ases 


Gas 

Ax 

Kr 

Xe 

Atomic ¥t . 

39.948 

83.7 

131.3 

l(eV) 

15.7 

13.996 

12.127 

I* ** (eV) 

11.5 

10.0 

8.4 

*a (cm^) 

1.64 X 10”^^ 

2.48 X lO"^^*" 4,02 X 10"“^^ 

Baa (dynes cm^) 

1,1 24x 10“^^ 

2.368X 10”^^ 4.0 X 10”"^^ 

Oaa (=“VeV) 

2.5 X 10“^° 

2.92 X 10“^° 3.75 x 10“^'' 

Cae 

7.0 X 10“^® 

8.2 X 10“^^ 1.0 X 10“"^® 


1.0 

1.003 

1.043 

**g+ 

5.6 

5.0 

4.54 


* Ref. 79. 

** Ref. 26. 



CHikPTER 3 


ElATHEIvIATIGil EOmroLlTIOE OE THE HROBIEM 


3 • 'I Boltzmann Equa ti o n : 


The Boltzmann equation is generally used for the descrip- 
tion of prohlems in kinetic theory of gases. This equation 
in'vol'ves the quantity f(x, c, t), the ^velocity distribution func- 
tion, which is the density of the particles in phase space (x, c,t) 
such that f(x, c, t) dx dc represents the number of particles 
at time t in the spatial ■volume dx centered at x whose velocities 
lie within the interval dc centered at c. The Boltzmann equation 


states that Bf/Dt, the rate of change of 'f’ along the free tra- 
;iectoiy of a particle is entirely the result of encoimters among 
the paiticles. For a group of identical particles this may be 
written in the form. 


Df _ 6f Sf ^ Z 6f 

^ a • eiE + Tr * 



coll. 


(3.1) 


where ^ is the external force acting on the particle, m the mass 
of the particle, and (6f/6t)^^^^ represents the change in f at 
a point (x, t) in space and time, due to collisions between 
particles evaluated on the assumption of dilute gases and 
binary collisions. 


The validity of the Boltzmann equation thus derived rather 
intuitively by Boltzmann (1896) has been reconsidered by many 
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authors like Y^von (1935), Bogolyuhov (1946), Kirkwood (1947)5 
Bom and Green (1949), Gross, Grad and ICrook (1960),^'^. As long as, 

(i) the gas is sufficiently rarefied, so that most of the 
collisions are binary, and 

(ii) the accelerations of particles over macroscopic and micro- 
scopic ranges are clearly separable, there seems to be no doubt 
about the validity of the Bolt2a:aann equation. In the case of a 
neutral and sufficiently rarefied gas, the low density and the 
short range forces between particles ensure the above two require- 
ments and the Boltzmann equation is valid. 

3*2 Validity of th e Bolt^zmann Equation for Ionized Gases ; 

In the case of ionized gases, however due to the long 
range coulomb force that exists between charged particles the 
above conditions are not satisfied v/hich casts doubts on the 
validity of the Boltzmann equation. The objection for the use 
of Boltzmann equation is that the assumption of binary collision 
does not seem plausible with long range forces and further the 
total cross section for collisions calculated on this assumption 
diverges. 

To study the justification of the Boltzmann equation for 
a plasma it is necessary to introduce some representative lengths. 
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Let j 

be the inter-particle distance, 
e2 

= _.™_ the landau distance, end 

V-T 

A = ( -i-™ ) te the Debye shielding length, 

^ STme"^ 

where n is the particle density, e the electron charge, k the 
Boltzmann Constant and 1 the temperature of the pla-sma. The , 
Landau distance gi-ves the distance bet?\feen two interacting 
particles, where the potential energy is of the same order a.s 
that of the mean themal energy of a particle and can be inter- 
preted as the mean distance of closest approach of the particles 
which is equivalent to the diameter of the particle. The Debye 
distance represents the effective range of two-particle correla- 
tions. 

The interactions of particles lying within the distance 
can be treated as strong interactions while particles lying 
outside Aj^ participate in weak interactions. Borces between 
particles separated by more than the Debye distance are statis- 
tically independent and can be taken care of by the average 
external force term in the Boltzmann equation. The forces arising 
from particles closer than A^j can be conveniently handled by 
the binary collision term in the conventional Boltzmann equa- 
tion. There is left the impoirtant range of distances between 
Ajj and Ajj (Pig. 3.1) in which the particles are correlated and 
yet they cannot be split into pairs for the purpose of binary 
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collision analysis. The number of such pa.rticles (— n4£) 

present within is lejrge, if n is large. ¥e introduce tb 

4-0 

dimensionless parameter g frequently used in literature as, 

5 -1 


Thus g is a representative of the ntunber of particles 
present in the Debye sphere. In most practical cases g is 

A iTN 

found to be small^ so that there are a large number of parti- 
cles within the Debye sphere and outside In this case, since 


M j3/2 


’^T. 5/2 


we have < Xj^ and hence weak interactions result from parti- 
cles lying in ^ < ^j). Thus the omitted range ^ 

is characterized by the super-position of many 'grazing’ deflec- 
tions. Treating these weak interactions as a stochastic process 
with many simultaneous Independent snail deflections, the result- 
ant contribution to the Boltzmann equation can be shown to be 
the Dokker-Planck tern which is a second order differential 

operator in addition to the integral operator which arises from 

21 

strong binary collisions . Making a 'grazing' collision appro- 
ximation, assuming a sequence in time of the independent random 
'binary' deflections, the Boltzmann collision term can be evaluated 

in the region ^ ^ ^335 yielding exactly the same fokker- 

21 

Planck term as in the earlier method 
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Ttie justification that has been given for the use of 
fokker-Planck equation is inconplete and a,ny analysis which 
demonstrates the validity of Pokker-Planck equation also seems 
to justify the validity of the hinary collision Boltzmann equa- 
tion in this round about way. Then the simplest procedure 
seems to be to include in the Boltzmann term ’grazing’ deflec- 
tions out to the distance instead of taking a sum of Boltz- 
mann and Bokker-Planck terms. Thus, though the validity of the 
Boltzmann equation for f-ully ionized gases is still open for 
investigation, the use of Xj^ as a cutoff distance in evaluat- 
ing the collision terms seems to be a, good approximation. (In 
the present investigation, g evaluated in the dovwistream equili- 
brium region takes values ranging from «-''0.01 to 1 and it may 
however take lower values in the earlier regions of the shock 
where the charge densities are much less than in the equilibrium 
region ), A simple method of obtaining convergent integrals is 
to use the IDebye-Huckel potential instead of the coulomb potential. 

5 • 5 Extension of t he Boltzm ann Equ ation for a Mixture : 

If the gas under consideration contains m different com- 
ponents, the Boltzmann equation can easily be extended as follows, 
under the same restrictions of validity. If f^(x, c, t) repre- 
sents the velocity distribution function for the j-th component 

(j = 1, 2, . . . m) then the Boltzmann equation for a, non-reacting 

1 5 

gas mixture can be written as , 
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ef . 

Tif 


+ C , 


ifi 

cx 


+ 


p. 

n . 

D 


ifi = 

0 c 


n 

S 

lc:=1 


J 


3k 


(3.2) 


■where is the collision integraJ. for scattering of particles of 

■typ© 3 due to collision with particles of type k. (k = 1,25 ...m). 
For a reactive riixture, the inelastic collision terms are to he 
included in addition to the elastic terns. The Boltzmann 
equation can he written for the reactive case aSj 


St 


4- C 


I- 

»i Tmwi i rrwO > nr,i» i n^ Mi 

OX m. 

^ j 


ef . 

e c 


n 

2 

k=1 




(3.3) 


where the superscripts e and i stand for elastic and inelastic 
collisions respectively. 


3*'^ Boltzman Equation for an Ionize d Gas; 

let fg_(x 5 c t), ^ 3 ( 5 ? £ "t) and fj^(x 5 c t) represent 

the velocity distribution fuQctions for atoms, electrons and 
ions respectively. Then the Boltzmann equations for the three 
species can he separately written as follows, 


Df 

= <f f >® + <f,f^>® + + <f f >* + <f f 

It aa ae ai a a ae 


+ <f^f.>i + <fife>'' 


(3.4) 


— = <fgf,>® + <4^ t 

^ a. 6X 8cl 6i 

+ <f f + <f f.>^ (3.5) 

aa ai 

Df. 

-TO- = 

+ + < Vi>i (3.6) 

where the significance of each term on the right hand side is 
discussed in the following. 
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5*5 Gollision j'e gas : 

In this section the significance of the va-rious elastic 
and inelastic collision terms is explained and integral forms 
of the -various terms ane given. A detailed evaluation of these 
integrals follows in Chapter 4. 

3.5.1 Elastic Collisions : 

In general, we write <f.f, representing the rate of 

J -K. 

change of f ^ , the distribution function of the particles denoted 
by subscripts j (at a fixed point x, along its trajectory) due 
to elastic enco-unters of j particles with particles k. k equal 
to j is also included, (jjk = a, e, i) . Schematically the process 
can be described in the following manner, 

h %'> * 

where c^^ and are the pre-collision velocities of the j-particles 
and k-particles respectively before the enco-unter and c! and cl. 

'"K I 

represent the corresponding velocities after the elastic inter- 
action. The symbol *A' represents the pa.rticle. 

A general expression for any elastic collision comprises 
of ' both loss and gain terms, because during enco-unters some 

i-particles disappear from the volume dc. centered at c. > 

’-^3 ^"3 

(in phase space) and some other particles of j type reenter 
this volume. Eor a process such as the one described above 
schematically in (3.7) the general collision term can be 
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written as 


'k1 


..(3.8) 

where is the elastic collision cross section for the j-k 

encounter. In general CT^]^ is a function of c^ , and the 

parameters that describe the encounter in the physical spaces 
■viz. the impact parameter h and the azimuthal angle G of the 
enco'unter; it also depends upon the interaction force law; 
g = is the relative speed of the pair of incoming 

particles k. for an elastic collision, as there is no 
exchange of energy between various degrees of freedom, the kine- 
tic energy of the particles in the cent er-of- mass coordinates 
remains invariant during a collision. consequently, the relative 
speed of the particles after the encounter is equal to their rela- 
tive speed before the encounter. Thus, 


S “ ^ 3-2 ".Skll " " -kl ! (3.9) 

Subscript 1 in is used for convenience. The variable with 

this subscript 1 is a dummy variable and disappears after the 
integration is perfoimed. 


3.5.2 Inelastic C o llision Terms ; 

(1) <f^f ; Represents the rate at which f decreases as a 

a, CL 3 , 

result of loss of atoms in an inelastic atom- 
atom collision of this type. When two atoms 
collide with sufficient energy, they may give rise 
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to ionization of one of tlie atoms resulting 
in the creation of an ion and an electron so that 
one atom is effectively lost. Schematically we 
can represent this type of collision as follows, 

MCa) + A'^Ccp + A(c^^^) + e(^^) 

where and ^c^^^ are the velocities of atoms 
before encounter and c'., c'>, c’. are the 
velocities after encoimter of the ion, atom and 
the electron respectively. The expression for 
the collision integral can he written as, 

^^a^a^ = - /^a^°a^^a^-albia^ ^ai 5 S 2. Soa 

i ... (3.10) 

where 0 - is the cross section relevant to atom- 

atom collisional ionization, and g = |c - c .I ; 

is the threshold value given hy = I. 

I is the ionization potential of the atom (more 

appropriately the excitation potentiaJ. to the 

^ m-imp 

first excited state of the atom) and ^ = - ■ - -t . 

- aa m-j+mp 

reduced mass of the system ( = "^ since 
m.j = 1112 = m^) . 

(2) <f f >^: Represents the rate at which f decreases result- 

^ 3 / 

ing from the loss of atoms as they get ionized 
hy electron impacts. When an electron homhards 
an atom with sufficient kinetic energy there is a 
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(3) 


possibility of the atom getting ionized resulting 
in an ion and an electron. The consequent dis- 
appearance of the atom is described by this loss 
term. This type of collision may be schematically 
represented by. 


where and are the pre collision velocities 
of the atom and electron respectively and c', c’^, 
Cg 2 3-^® 7he po,st collision velocities of the 
resulting particles viz. ion and the two electrons 
respectively. One electron plays a sort of dummy 
role in the process similar to that of an atom in 
<fafa>^. expression for the collision integral 

describing this process can be written as, 

=-Aa(2e)^e(5ei)a5,e « « > 8oe 

where is the collision cross section for the 
ae 

atom-electron collision resulting in ionization 


e1 


g^^ is the threshold 

06 




of the atom and g = | 

TOlue given ty ^ n g^e = ^ae = ( 

being the reduced mass of the system; 
since << m„; ni and m^ mentioned above are 

6 cl 8. 6 

the masses of atom and electron respectively. 


Represents the rate at which f^ changes due to 
collisions between atoms and ions resulting in 
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the ionization of the aton. Schematically this 
process can he represented hjj 

- A+(cp + A*(o!^) + 3(0;) 

where c and c.. are the precollision velocities 
of the atom an 5 ion respectively and _c^5 are 
the post-collision velocities of the nev/ly horn 
ion and electron respectively and c!^^ is the post 
collision velocity of the dniniriy ion which parti- 
cipated in the collision. This term has been 
neglected in the present studies; (see Section 4.. 4.1). 

(4) <f. f„>^; This represents the rate at which f increases due 
to the recombination collisions between ions and 
electrons, resulting in the formation of new atoms. 
This process is the inverse process of ( 2 ) above in 
which three bodies participate. A third durray 
particle, an electron in most cases, serves the 
purpose oi transferring the required energy in the 
process. Schematically the process may be repre- 
sented as follows, 

A'^(c|) + e^(c^) + 02(0^1) - A(c^) + ©2(0^^) 
where c„ and c^. are the velocities after colli- 
Sion of the resulting atom and the dnmmy electron 
respectively; c^, 


are the velocities before 
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collision of the incoming ion, the electron to 
he captured and the dummy electron respectively. 

Only those conihinations of _c|, _c^ and that 
result in c^ for the final atom are considered. 

The contributions arising from the recombination 
processes can be obtained from the corresponding 
ionization processes (see Section 2,6). 

(5) <fQfg^>^ : Represents the rate at which f^ changes due to 

collisions between a,tons and electrons resulting 
in the ionization of the atom. In this case there 
are loss as well as gain teims. If the electrons 
having velocities in around participa,te 
in this type of collision, they villl be scattered 
out of their original phase volume causing in 
effect a loss term for fg. On the other hand, I 

& few of the electrons crea^ted in this process ' 

may have velocities lying in dc^ around c„ appear- 
ing as a gain tern for f^. The process may be 
schemat ica,lly represented as follows, 

I 

a) gain: + ei(Cg^) + egC 

b) loss: A(o^,) + e^(Cg) - A^(c!,) + e^(o^) + 82(04 

In (a)_c^^ and c^^ are the velocities of atom and | 
electron before collision v/hich result in the crea- i 
tion of a new electron whose velocity lies in 
around and^^^ are the velocities after ! 
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the collision of the ion and the dinnniy electron 
respectively. In (h) and are the velocities 

of the atom and the electron he fore collision and 
c ’.^5 c* and c'. are the velocities of the ion, 
the dniany electron and the newly horn electron, 
respectively, s,fter the collision. The collision 
integrals may he written separately for gain and 
loss cases as follows, 

"^^e^a^'gain ' /^a^«£a1 ^ ^e^<£e1 -b^e ®'^eV 

S’ > Soe (3.12a) 

"^^e^a^ 'loss j ae^ ‘iSal ’ 

g > fioe (3.12h) 

where cT^g is the cross-section for atom-electron 
collisional ionization and g*, g the relative speed 
of the incoming particles. 

(6) <fgfj^>^ ; Represents the rate of change of f^ due to elec- 
trons participating in the ion-elcctron recomhi- 
nation process. This contains hoth gain and loss 
terms. The loss tern appears because an elec- 
tron disappears in each such collision. If the 
dummy electron which participates in the collision 
has a velocity before collision lying in 
around c„, then this can also he counted for the 
loss tern. If the dummy electron acquires a 
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velocity lying in at alter tlie collision 
it is to te treated as a ga,in term. The process 
is s ch em at i easily represented as folloT/s; 


(a) gain: A^( c[^ )+e^ ( c^^ )+e 2 ( ) 

(b) loss: A'^(Cj_^)+e^(^)+e 2 (cQ^) 

4^(Ci^)-t-e^(Cei)+e2(Ce) 




wTaere in (a) ^j_-|s ^@1 >£4 velocities 

before collision of the ion, the electron to be 
captured and the duaiiny electron respectively; Cg_-] 
and^e velocities after collision of the 

atom and the dummy electron respectively; in (b ) 
c -^5 c^ and c„. are the velocities before colli- 
Sion of the ion, the electron to be captured and 
the dummy electron respectively where as and 
c^^ are the velocities after the collision of the 
atom and the dummy electron; in the second equa- 
tion of (b), > Xel ^^'^vXe veloci- 

ties before collision of the ion , the electron to 
be captured and the dummy electron respectively; 
^a1 ’ S.e ^®P^®s®n't ■fche post collision velocities 
of the atom and the dummy .electron respectively. 


(7) <f^f^>^: Represents the rate at which f^ increases follow- 

ing creation of new electrons during an at on- atom 
collision resulting in ionization. The process can 



"be schematic ?lly represented as follows, 


A(op+A(c4P - a" 


(£i)+A(o^P+e(Og) 


where and _c^^ are the velocities before colli- 
sion of the two incoming atoms; Se 

the velocities after collision of the ion, the 
dninmy atom and the newly hom electron respec- 
tively. Ihe expression for the collision integral 
can he written as 


<^afa> 


/^a^^a^^a^-al ^CTaa ^'^a1 


g‘>Soa (^•■'5) 

where (T^ is the cross section for atom-atom 
collision that results in ionization and 
g' = 

Represents the rate at which f changes due to 

w 

generation of electrons followed hy atom-ion 
collision resilLting in ionization. This term is 
ignored following the arguments given in Sec- 
tion (4.4.1). 


Represents the rate of increase of f^^ due to ions 
being created in the atom-electron collisional 
ionization process. Schema.tically this process 
can be represented as follows, 

A(=a)+ei(S^l) - A'"(Oj^)+e^(CgP+e2(Og2) 

where c’ and c'. are the velocities of the atom 



and the electron lefore the collision and c. , 

-Se1’^e2 velocities of the ions "tlie 

two electrons after the impact respectively. The 
expression for the collision integral may he 
written as 




1 ae ® * ^e 1 


g’>-goe (5.14) 

where CT ^e cross section defined in (2) 

above s and g' = |c^ - c'^^l . 

Represents the rate at Tfhich f^ decreases due to 
the disappearance of ions in the ion-electron 
recombination process resulting in a neutral atom. 
This is the three body interaction process des- 
cribed in (4). Schematically the process can be 
represented as. 


f (ai)+®(=e1^+®(£e2) ■* 

where and _Cg 2 are the velocities before 

collision of the ion and the two electrons; c’^ 
and are the velocities of the newly formed 
atom and the dummy electron respectively. 


Represents the rate of increase of f^^ due to the 
creation of ions in atom-atom collisional ioniza- 
tion. Schematically this process may be represented hy 
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(12)<faV'^ = 


where and are the velocities before impact 
of the two incoming atoms; c . . are the 

velocities atter collision of the ion, the d-ummy 
atom and the electron created respectively. The 
collision integral may he written as follows, 


<V^>' 

8/ Si 


■/ aa ®'^a1 


8' > Soa (3-''5) 

where CT ^ is the cross section defined in (1) 
and g' = lo4-o;,|. 


Represents the rate at which f^ changes due to 
any possible creation of ions in the atom-ion 
collisional ionization process, following the 
arguments of Section ( 4 . 4 . 1 ) this tenn is neglected. 



CmiP'TLR 4 


GOVERra& DIFEEEIiHTIJOl' EQUATIONS 


In this chapter suite^ble functional forms for the distri- 
bution functions referred to in Chapter 3 for the three species 
and the appropriate moment functions are chosen. For the choice 
of these and the collision cross sections frorxi Chapter 2, the 
collision integrals of the integro-differential moment equations 
are evaluated analytically so that the system of equations 
reduces to a set of nonlinear first order ordinaiy differential 
equations. The boundary conditions relevant to the problem are 
also discussed, 

4 . 1 Distribution F -u nctio n _f or _ the G^a, s : 

It is well known that a bimodal distribution function of 
Mott-Sniith type con adequately describe the nonequilibriun 
region inside a strong shock v/ave even though it does not re- 
present an exact solution of the Boltzmann equation even at 
infinite Mach nuinber^^. Apart from its simplicitjr there are 
other reasons to continue the bimodal ansatz for the description 

of the shockwave. The exact solutions of the B-G-K model equa- 
1 7 

tion suggest that the true distribution function does resemble 
the Mott-Smith’s ansatz although later detailed asymptotic 
studies show that it is unlikely to be an exact solution. 
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Inanthasayanam and Narasimha^ have' shown for the B-G-K model 
that the hinodal ansatz can gf-ve reliable estimates of the 
density-slope shock thickness presided suitable criteria are 
employed. Recently Bird employed a direct computer simula- 
tion Monte-Carlo technique to study the stmeture of strong 
shock waves. The velo.city distribution function within the 
shock wave has been illustrated by computer display photographs 
with the molecules represented as dots in the velocity space. 

This result of Bird gives qualitative support to the biniodal 
model. However, Bird's results predict a certain .degree of 
aS3niiimetry in the density profiles contrary to the bimodal predic- 
tion, though the ansatz can predict reasonably the maximum- 
density shock thickness profiles. 

The bimodal ansatz is more general than that originally- 

q -i 

given by Mott-Smith in as much as the parameters describing 
the two Maxwellians are functions of x. It can be shown that 
in an inert gas problem, in which the collision terms of the 
differential forms of the equation of continuity, momentum and 
energy are zero, the choice of these parameters is not comple- 
tely arbitrary I and if the parameters of one of the Maxwellians 
are considered as constants, then it is imperative that the 
other two also should be constants for consistency. However, 
when the inelastic collisions present make the collision terms 
of the conservation equations nonzero, it will be possible to 
consider the parameters of one of the Maxwellians as constants 
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corresponding to the upstrean equilibrimi conditions, keeping 
the parameters of the other Max\'''0lli8jis as functions of x. 

The original physical idea of Mott-Sniith namely, that 
the state of a point within the shock v/ave can be described 
by considering the point to be characterised by a certain 
number of molecules of upstream equilibrium region and a certain 
number of the downstream equilibrium region has been extended 
by Haight and Lundgren^'^"^^ to the study of shock ^eye structure 
in diatomic gases taking the excitation of rotational degrees 
of freedom into account. Haight and Lundgren considered a 
trimodal distribution function, two Maxw^ellians of the three 
being characterised by one upstream and one downstream equili- 
brium regions; the third Max^/ellian is characterised by an 
intermediate equilibrium state defined by permitting only the 
excitation of the translational degrees of freedom, while the 
energy in the rotational degrees is unchanged. The trimodal 
'distribution function is then a sum of throe Maxwellians with 
weighting functions dependent on x, and the parameters characteris- 
ing the Maxwellians being constants. The hypothesis of a trimo- 
dal distribution function is essentially an extension of Mott- 
Smith’s bimodal theory for a monatomic gas and consequently is 
most applicable to strong shock waves., - far the trimodal 
ansatz has only been sporadically accepted. In what follows, 
a trimodal distribution function for the atom gas is assumed 
where the two Maxwellians are characterised by the two equilibrium^ 
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regions upstream a,nd downstream, and the third is characterised 
hy an inert gas equilibrium state where only translational 
degrees would he excited. The atom distribution function can 
be ?/ritten as, 

ij. = E n.(x) f. I 0 = 1, 2, 3. (4.1) 

- ;3 3 3 

where n^(x)'s are the vi/eight functions, and 

V 

exp [-(Cg_-U.i)VPjl ^ (4-. 2 ) 

is the molecular velocity for the a-toms, and Pj is defined by, 

p. = 2kT^/m^ . (4.3) 

In eg., (4.3), is the mass of an atom and k the Boltzmann 
constant, and T- in the above expressions represent the mass 
velocities and temperatures of the atom gas. 

Subscripts 1 represent the quantities corresponding to upstream 
equilibriuiii condition, 

2 correspond to intermc-diatc fictitious state v^^here 
only the translational degrees get completely 
excited (corresponding to the downstream equili- 
brium in an inert moncatomic gas) and 

3 correspond to the downstreoj;' equilibrium state; 

"5^ is a unit vector in x-direction, representing the direction 
of the flow. 
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4.2 Ele ctron cud Ion functions; 

Petschek and Byron° have e^rgned that a tiaxwellian 
distrihution can be a-ssumed to exist in an electron gas if 
an electron ha.s appreciably r-ore tha,n one collision with another 

electron for eveiy 10 collisions with etons; this requires 

—8 

0 . degree of ionization grea/cer than 10 . Within the realn of 

the present work the electron gas as well a.s the ion gas afe 
each assumed to ha.ve Maxwellian distribution functions, v^rhich 
can be written as, 

^e " ^e^^^ e^cp [“(Cg-Ug (x)i) Vpel 

and f^ = n^(x) (tcP^(x) exp [-( (x)i) J , ( 4 . 5 ) 

where c_ and c- are the absolute velocity vectors of the elec- 
trons and ions respectively; p is defined by 

Pp = 2kTp/mp , p = e, i . ( 4 . 6 ) 

Subscripts e and i stand for electron and ion respectively; 

n (x) 5 U (x) and T' (x) are the number density, ’mass velocity 

P P P 

and temperature of the p species respectively and 111 ^ is the 
mass of the p pa.rticlei 

^ • 5 Moment Pun c t io n s ; 

It has seldom been possible to solve the Boltzmann equa- 
tion for the distrihution function f, and most of'ten only a, 
set of moment equations obtained by talcing moments with suitable 
functions 4) (c), of the Boltzmann equation have been solved to 
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describe the irnknown parasieters in the assumed fom of distri- 
bution function. Amongst the choice of moment functions 4>(c^)s 
(})(_c) = 1, c^ and c result in the fai:iiliar con servant ion equations 
of na.ss, moment imi and energy. The difficulty that ha„s been faxed 
todate right from Mott- Smith’s vrork is the insufficiency of the 
number of equations to determine all the unknoras of the problem 
with these moments. Mott-Sliith used an additional 4> (^) to 
complete the solution by choosing c^ and c^-^ in his analysis. 
Attempts by Sakurai"^”*, Rosen^^ and others did not yield much 
so far as eliminating this ambiguity is concemedjas pointed out 

by G-ustafson and Narasimha/ , Neteitheless it is felt that 

2 

the choice of c^ seems to give closer estimates of the density 

profile shock thidkne sses, though there is no rigorous justifica- 

2 35 

tion for the choice of c„ . Marasimha and Deshpande frora 

, Jx. 

their minimum error solutions have recently concluded that 

2 

results should be better than c^ particularly for higher Mach 
nuDibers, However, their method is not free from certain arbi- 
trariness either, since the definiti.-'U of the 'error* itself is 

2 

not unique. The present work continues with the choice of c^ , 

from the point of view of simplicity over the minim-um error methods. 

15 2 

Bird's results also seem to support to the choice of c^- . It 

should be noticed, as highlighted by Rode and Tanenbaum that 

the solution is highly sensitive to the choice of <{) (c). 
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^ G-ove ming Diffe rential Eq- uations ; 

Expressions (4.1), (4.4) and (4.5) can be substituted 
in the equations (3.4), (3.5) and (3.6). Using the linear 
cross sections froii Sections (2.5), (2.6) and (2,7) -various 
moment equations can he derived. 


Ui miher Oon servat ion E qimt ipn_ for Atpn s j 

As the elastic collisions do not result in any change 
in the number of particles of any species, the elastic colli- 
sion contribution is zero. The atom-ion collisional ionization 
is neglected through out, tho-ugh, this cross- section is of the 
same order as atom-atom collisional ionization cross- section. 

This is justifiable because in the atom-electron collision 
dominated region it is as insignificant as atom-atom collisional 
process and in the atom-atom collision dominated region the ion 
population is too far less to compete with the atom-atom process. 


Integrating equation (3.4) over the c space, we get, 


4 


‘dx [ J ^j ^j ] k 

where the are as given belo-w. 


(4.7) 


A(j)-| = rate at which atoms are lost due to atom-atom 

inelastic collisions, 

- /fg_fa-[G"aa ^ ^a iSa1 ’ ^ ~ ^oa * 

The factor 1/2 is included to account for the duplication of 
collisions. Using, 
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^ ^aa ■^os^ “ 1* '^aa , (4-* 8) 


^'^1 " H ^-a °aa/|^ ^D j [“( ^/P jl] 

•[2 exp [-(c.r-Uki)Vpk]] (g^-gia)g 5Sa^^ 

Using the results of Appendix A- I, 

1 


A(|), = 


I'V C, 
lefTt ^ 


^iV^Pr ^1 ( U. 9 ) 


where , 


:Pi(U,Y) = ^fn [srf(Y+P)-Erf(V_U)j + (| » ^ 

+ exp [-(U^+V^)J (2U^+5) cosh 2UV 


+ (2TJ^+3) 2 sinh 2TJY 


(4.10) 


and, 


i'-](0, Y) - 4(2 + Y^) exp (-Y^) ; 

and p.^ = p. + 


(4.11) 


4^2 rate at which atoiiis are lost due to aton-electron 
inelastic collisions. 


IVeC^Q^e C £ 11 Sq© 


^ ae 
A(p 2 


.2 „2 


2-1^ ae "ae ^^""§ 06 ^ “ 2' "^e ^ae ^S^-^oe), 

?“e exp [-(o^.u.t) tp^] 


• (g ~&oq) g • 


Using 
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Using the results iron Appendix A-I for the integrals, 
the expression for A ^2 witten as, 


whe re , 


A (j) 


Ue 


1 


Dg C n z Fi C-^-- , ) 

Wtd S’ a. e 3 UZ'- ' ^rp.^ Tfp. 


u. 


je 




= P3 + Pe 


(4.14) 


A 4) 


A 4 


rate at which atons a,re created as a, result of 
re comb in at ion between a.ton-ion~ele ctron, and 
can written froin ects. (2,34) as, 


- A(i) 


1 


4tc 


3/2 £a ,3/2 


n .n - 
1 e 

nr"* 




a 


(4.15) 


rate at which atoms are created as a result of 
ion-electron-electron recombination, fron 
eq. (2.34), 


~ A4) 


ST „'372- (8oe/Pe) 

e Pg 




n 


a 


(4. 16) 


In the eqs^ (4.15) and (4.16) 


2kl 


n. 


a 


E n.; 

D ^ 


a 




a 


n 


a. 


rn 


a 


and is the atom temperature defined in eq* ( 5.21 


. 13 ) 


Substituting expressions (4.8), ( 4 ,l 3 )g (4.15) and (4*16) in 
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eq. (4.7) s the a^toi-i conservation equation can he written as 


(2 n Tj ) = „ . 

ax D D 


•] 


n 0 

“a aa 


hte3/2 h .3/2 1 exp (22?-..) ^ 

bu- -T ^ 2|3 


+ 


7^ 




oa 


p^/2 j, ( _aiL , 


) 


1 

4'ra 


■'e *^ae 


b4.^/2 & « 


1 


/^oe^ ^i-^e 

p- v’ ™'5T2 ®^(~n"-) 


U-, 


2 n. P, ( .iLe, ^ 

3 ' ' '''ho 


"oe 


'he 


) (4.17) 


4# 4# 2 Cv- Moment Equation for the Atom Gas: 

) a . 


? w (n? + i Pi) 


2 

6 

Z A 0. 


Multiplying eq. (3«4) hy and integrating^ we get ^ 

q, P (4.18) 


k=1 


whereAGj^'s are described belov\f. 

'' ®1 =7 /=L (k - Va4 e + 

where b is the impaxt parameter a.nd 8 the azimuthal angle. 

Using the results from reference 18, (4.19) can be vyritten as, 

A 9.. = 4 r (c’t - G^^) f„ g bdb d£ dc dc.., (4.20) 

1 2J ^ ax aX'' a a1 ® -.^a'-^ai 

Substituting for f fron eq, (4*1), and using th^ results from 

a 

Appendix A-V, eq. (4.20) reduces to. 


.19) 
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A 9 


B 


■ ^ 2 ( 5 ) E E n.n, 

2Y2 °a ^ 0 k 3 ^ 


V 


( 4 . 21 ) 


where B „ is the force law constsnt for the ctoia-atom interaction. 
a,a 

''02 =/4 (f; - fa fe) S -fS fca 4?e 

= /(4 - 45 fa fe S I’" '•'= fca S=e 
from the results of reference 18. 


Inserting the esipressions from (*1.1) and (4.4) for f^^ 
and fg, A 02 can he v/ritten as 

^^02 = / h exp [-(Og^-D^5)‘'/Pj]|ng()tpg)'^/^ 

• exp [-(c^-U,l)4pe] (<=^ 1 - 4 ) 8 * ISa 

Evaluating the integrals using the results of Appendix (A-Y) , 

A 02 can he shov/n to he equal to, 

Ae2 = e n S n p P (^^- , -il • , , 6) 

(4.24) 

where , 

P ^ ( U, Y, X, V) = ( ( 5 ) ~ I 2 ' 5 ) ) " 

- 2A^(5) UX+A^(5) (W-Y); (4.25) 

A^(5) and A2(5) are define... in Appendix. 

* 83 = /4 fi - A fi> s de dp^ do^ 

= i <4 - 4 ^ t. A 8 .db dp^ dp.. 


(4.26) 
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Using the expressions iron (4.1) and (4.5) for and 
fg, A 9^ can he written as. 

Ad ^ ““5/2 * -v~v. P TT *l1 nn { d 


A 0^ = ( 
0 


S n^(ixpp exp [-(Cg-U4)Vp;5l j n^(xpi)' 

D 


exp [~(c.-uni)^/i3j (c’|-4j g hdh di (4.27) 


which, after using the integrals f roi:i Appendix (1-V) reduces to 
the forn, 

„ .1/? - _ n, p, U., 1 


fZu (■^+- 

ct 


B • S pj: * 7 ^ ( 5 9 ^ ^ 

i ^ ^ Pni 'rPii " 


31 ^31 ^31 

(4.28) 


In obtaining the elastic collisional expressions (4.21), 
(4.24) and (4.28) the inverse-fifth-power laws of force described 
in Section (2,7) have been used, 

A = - ■g' fg_.| (j g^Q_ g ^a. '^a1 ’ S > goa (4.29) 

Using the distribution function described in (4.1) s^dCT^g^^ 
from ( 4 . 9)5 2q. ( 4 . 29 ) can be written in the forin. 


B Caa/l:2 ^^3 ^^^3 exp [-(c^^-u^f) 2/p . 




'""a 


1 w jsysALaiwi * r s» hl x« | 

‘ 9 \ ? Q 9 9 / * 

*" aPn aP/^ iPi^ 


(4.30) 



.,5/2 

li . Iii„ [ 

3 k t 



,Pf 

in,. , 

. Z'i 2-„ 


'^Pjk 


(4.31) 
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after using the results fron appendix (A-Il), where, 
I'2(U5V,X,W) = ][% |^Erf(¥+X)-Erf(T/-Z)j V+U^)X^ + | UVX" 

3 ttO 1 / rr m2 ^ / 1 tt , tt2 ^ / 


+ V 
1 


+ ^-(3-W^)(^ V+lT^)j X + 2 2nr(3-Yf) + 


+ ‘gr V+t7‘^)(3-2W^) IjT(3-2W^) 


1 IT? / -7 n .7r2 ^ 1 


1 ^ 


TO- 2 \r) 


,-3 


2^, 1 
2,..2. 


+ exp [-(Yr+X^^)] 


• j^ J(2' Y+U^)X^ + 3UVX + (^- V+U^) + I 


+ X UY (3+41<’f^) “ - Y^ (3+4W^) 4. 

X • X"^ 


cosh 2 IX 


+ |^('gY+U“)X^ + 3UVX + -2’ V^(3+2W^) + | (^Y+U^) 

- I UY ^ -4 I I sinh 2wl] 

d X 'r j ^ j 

and, 

E2(U5 0, W) ” .16U^ + 8U^1Y^ + + -| + 6 exp (-W^) 

(4.32) 

* ®5 ' - /^a <T^e S JSa i=e = S > g^e - 

(4.33) 

Substituting for f and f from eqs. (4.1) and (4.4) and using 


0* defined earlier, 

Clt/ 


^ 6, 


■ 2 ’ -b [t [-(Cg,-U^^)Vp 3 lj 


• ^^(TcPg) exp [-.(Og-U^i)VPe] 

. (g^-goe) S 4^ 


-^n2 


r . 


" “4 Poo “e ? ^0 Pjo^ 

D 

U. p. U. 

. E.-. f iL. .-^9?. 

^je ’ Pje ’ ^Pje ’ ’^Pje 


be 


) 


(4.34) 
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after evaluating the integrals using the results from 
Appendix (l-Il) . 


A e 


6 


of„ 

6t ^rec. ax 


.2 . Px. 


which can he shown to he appro xhoalely equal to (U^ + • 2 ~) times 
the recomhination rate. Therefore, 

3 


A Bg = 471; 


°a «' 


=+ n 


ri^.n^ (Uj^ + -^)(2 


\j n .11 

z ae 1 e 


.2 

=‘oe 

Te 


■) (4.35) 


The contribution from atom-electron- ion recombination is neg- 
lected. 


Substituting for A0j^’s, eq.. (4.18) reduces to, 

S-.V2 ^ ^ 2 

A2(5) Z 2 n. n. 
a j k 3 


Z n.a.(u2 +1 p.) 


TC 


/ aa\ 


®ae V2 

.,2. (-^) en^Zn.p.^P, (^^p. 


je ^de ^e’ 


+ iTSitC--) n. S n. ^ ) 


u. 


P-1 u,. 




"oa 


S'lit 

1 




7r,i ’^e C,,e “e ^ M £^> ifir~) 


u. 




g. 


oe 


+ 4n 


"a h 


Pl .. 


ri 


-2 °ae - 2 -) (2 + ) 


=oe 


'e 


'3® *"^30 

(4.:;6) 


4.4.3 Humber Co nserv ation Equatio n for kl^ct :ron_s_: 

.. Integrating eq (3.5) over the c sp8.ce, gives, 


dx 


(Ve) = 


4 

2 A 41^ 


k=1 


(4.37) 
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v/here the right hand side of eq. (4.57) is the negative of that 
of eq. (4.7), since for each aton lost or gained there exists 
an electron gained or lost irrespective of what the velocities are. 


4.4.4 Momen tini Equ ation for t he El ectron Gas ; 


The electron nocentun equation can be got by nultiplying 
equation (3.5) by c„^ and integrating over c^. The left hand 
side of eq. (3.5) for x-direction for s'fceo-dy ca,se is. 


ex 


efg Gfg 


(+.38) 


where is the product of the electric field in x-direction 
and electron charge = Jlie. Multiplying (4.38) by c and 

"BIX 

integrating over c^ yields that the left hand side of the 
moinentuni equation 



since (f„)^ = 0. 

^ e oo 


Moment UC .1 equation is given by, 



The electron-electron elastic interaction can be assuned 
to result in no net change of moinent -uri . 

^a "e § 5^e ^a‘ 

Using the results from Appendix (A-III), after substituting for 
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and fg from (4.1) and (4.4), eq. (4.40) can be written as, 

— \ 

= 2 % A-,(5) n^ 2 n^ U.g ^ (4.41) 

e 2 

*'*'2 = / (“ex-^ex) dCg (4.42) 

wbicb can be sbovm -using the results from Appendix (A- III) 
after substituting for f^ end f^, from eqs. (4.4) and (4.5) 
to be equal to. 


AiJ)2 - -^1(2) n^ng 

e ''Hot 


'01 ^ei 

exp [-(UQj_/pgj_)]- Erf (:;-r-|— -) 


ei 


ei. 


(4.43) 


whe re . 
Appendix. 


= Pg + "'^Gi “ ^"^6“^! A.j( 2) is defined in 


Each electron produced or lost in inelastic collisions 
can be assumed to hruwe an average local mean speed Ug and a 
random speed corresponding to the local electron temperature 
I‘ . Undi^r this assiiraption the contributijn to nonentuni by 
inelastic collisions, per electron is -unaltered and hence these 
terms can be omitted. A detailed quantitative estimate of 
these terms made but not included here justifies the assumption 
of negligible contribution to the momentum equation from inelastic 
collisions. 


Hence the electron momentum equation can be written as, 
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4- Cne(Ue' - ?)] 


+ 271 % A^(5) S U.^ 


4f7i; (-^r A. (2) n. n 

^-‘e ' ^ ^ ''Hie ^ei 


exp(-(^«"0) 

Pei 


Si)) - jTi!. 'f^ol (^ei.) 


w;iii 

(4.44) 


4.4.5 E lectron Energy Eq uation : 

Tlie electron energy equation can "be obtained by talcing 

p 

Cp -moment of eq.uation (3.5) which can be written as 


dr- , 2 ^ -r 2n U Ee 3 

Pvn TT f 11^ 2. n M ^ ® ^ — V 


He lie i Pe)] + 


Axk ■ ('’••+5) 


The elastic collisions between electrons do not contribute to 
the net change of energy of the electron gas. 


Ax 1 


= f(c2”cf) f„f_ g bdb d£ dc^ dc_. 


(4.46) 


which, after using f and f^ from (4.1) and (4.4) and using the 
results from Appendix (a-IT) reduces to, 

®ae VP Up P. U,. 

AX 1 - 2” --iC5) "e ^ Uj Pje P^’ 


where 


Pp (UjVjX,!) = 3(W-Y) + 2W1^ - 2UZ 


(c'^ - c^) f. g bdb d£ dc^ dc. 
'• e S'' e 1 ° ^e 


(4.47) 


(4.48) 

(4.49) 


which reduces, after using the results of Appendix (A-IV) to 
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Ax 


471 A^(2) n^iig 


P (>_ is. 

5 Wei’ Pel’ 


u, 


ei 


TFei 


, 1 ~£) 


(4.50) 


where, 

P 3 (U,V,X,¥) 1 Erf(X)(W _ |) + exp (-X^) (J - Y) 

(4.51) 

With the same argument as gi'ven in deriving the electron 
momentum equation, it can be assumed tha.t the change in electron 
energy due to processes (5) and (6) of Chapter 3 is mainly due to 
the energy I, being lost from the electron during ionization or 
being gained by the electron in recoinbina,tion, and hence it can 
be written that, 

2 

Ax 3 = gj (Recombination ra.te - Ionization rate) (4.52) 

2 

where, g^ = 2l/n^. 


I'hus the electron energy equation can be written in the 


form, 


^ DC 2n^R Se 

^ TT Z' I 5 N“1 GO 


m. 


B ^ 1/2 
• + 111 (-»■■) A-|(5) n^ 




-Tx C^e^e^^t + f Pe)] 

• ? Y Pje Hj’ r// 

H- 4„(^)2 4,(2) P3(^% ^ 


Yl 


1-e) 


ei 


1 


w Pe 


n . n 
i es 




. S n 


.3/ 


0 "^3© 


1^ fpje’ ^ 


3© 


(4.53) 
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4.4.6 ConsejCYsjb i£n Equations for -the Plasma : 

Conserva.t ion equations for tlie plasma 8,re obteaned by 
the a^ddition of the conser'va.tion equations for the three species. 
The overall conservation equations of nassj moment uai an,d energy 
v/ill be almost exact since the approxiiration .errors in the colli- 
sion integrals of the component equations arc .nullified by the 
addition of the component equations. 

Multiplying the mmiber conservation equations for o.toms, 
electrons and ions by their respective masses and adding yields, 

■h- L “a ? U. 4. B. n. + ng Uj. ] = 0 (4.54) 

which represents the conservation of na.ss for the plasma. 


The nomenturn equation for the plasma can be obtained by 
multiplying eqs. (3.4), (3.5) and (3.6) by c^^, nig c^^ and 
m. c.„ respectively, integrating over c^ and c. and adding. 
During a collision, elastic or inelastic the total monentum of 
the participant particles is conserved in the absence of any 
extenial force. The only external force present, namely, 
electric field arising from possible polarization has already 
been included in the force tern of the Boltzman equation. The 
resultant change in momentun of the entire gas due to collisions 
is therefore zero so that the momentun equation is, 



+ Ee(n--n . ) 

“ 1 


(4.56) 


0 
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where the last teEJ in (4.56) arises iron charge separation. 

The energy equa,tion can he obtained by e-va,luating 
12 12 12 

? ^^e°e 7 . ”i°i noHie-nts respectively of the eqLia- 

tions ( 3 . 4)5 ( 3 . 5 ) and (3.6) and adding. The elastic collision 
terr.Ts disappear, since the tot .--I energy is conserved during 
every elastic collision. In the inelastic collision an amount 
of energy equal to I is effectively lost or gained according as 
the process results in ionization or recoinbination, and the rate 
at which energy is gained can be written a.s 


AE = I (Recombination rate - Ionization rate) (4.57) 

Substituting in the right hand side of eq. (4.57) from 
eqs, (4.37) and (4.7) for the net rate oi production of charge, 
eg.. (4.57) reduces to, 

''S = - I 4 = - •? =>8 Si ci (mOo)- 


Using this, the energy equation for the plasma can be written as 

4: “a f + i Pp + S ! Pe> 


+ ^ m.n.U.Ca^ + I pp] + eE(n^Uj,-n.np+ ^ g 


■ 4c (•'Wo) 


0 . 


(4.58) 

Subtracting the electron conservation equation from ion con- 
servation equation yields, 



(4.59) 


since the rates of production of ions and electrons are sarie. 
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Intef^rat ing (4.59) and using tTae condition fhat, 

X - - oo, = n. = 0, yields, 

-L 

° ( 4 . 60 ) 

Equation (4.60) also inralios that the net current flux is 

zero at eveiy point when there is no externally applied electric 

field. 

Using eq. (4.60), eq.(4.58) con he written as, 

"dx l9*^a ? °e^e'^e^"^e 2' ^e 

0 

+ m^n^u^(uj + I P^)] = 0, (4.61) 

from which it is seen that the electric field terri does not affect 

the energy equation. This is because in the absence of external 
fields, the sum of kinetic and potential energies remains unaltered. 

G-auss Equ ation; 

lastly, the electric field E is described by the Gauss 
3 quation, 

V. E = 47ie (n.-n ) 

i^hich for a one -dimensional problem reduces to, 

^ = 47te (n^-Ug) (4.62) 

• 5 Nond imensionalization o f the E quations ; 

The various physical quantities of the problem can be 
ondimensionalised by using the upstream n^ and as follows, 
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% ’ 


U. = '3. = fj. -i 


*1^2^ _^t)0jXa 


X = 


\/ ^ -| 
X 


E 

W" ^ "" L 


(4.63) 


a 


where L is a suitably chosen length scale. 


I'req.uently the nic^n free path x , corresponding to the 
oolcl side of the gas hs.s heen used for nondinensionalising the 
shock thickness and the equations. This neon free path, unfor- 
tunately is "vciy sensitive to the interro locular potential law 
shosen in the kinetic theory nodel or equivalently to the visco- 
sity-temperature law chosen in the continuun model. For values 
)f s, the power in the force law P = B/r^, > 1 , X for 

L| ->• 0° and hence makes the nondimensioncl. shock thickness go to 
;ero though the actual thickness of the shock remains finite. 

'rom time to tine, the need for a more nppropriate length scale 
nsensitive to the choice of the force la,w has heen stressed, 
ighthill''-^ has recognized in the solution of ilavier-Stokes 
quations for the shock wave thickness, that the downstream 

ean free path does not go to infinity as il-| ^ 0° for any vaBue 

42 

f s. leipmann et.al. have shown that the use of the mean 

ree-path at the sonic point inside the shock is particularly 

52 

ppropriate for the Havier-Stokes solutions. Mukenfuss from a 
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study of Mott-Snith type solutions suggested fhe nean free path 
corresponding to a point inside the shock, where v = -g' (Ref. 51, 

= n^/n^). The l 2 .tter lengths, where a knowledge of the con- 
ditions at an interior point in the shock is required are diffi- 
cult to use in practice, IxTcrasimha. and feshpande^^ suggest the 
mean free-path side of the gas in their study on 

Eiinimun error solutions of the Boltzmann equefion for shockwave 
structure. Xg always rema-ins finite Both in the continuum and 

kinetic theories, further, they shov/ that the temperature 

-1 

varia,tion inside the shock after the point where v = -g- is not 

significant for large Mach numbers so that the mean free -path 

on the hot side of the shock is of the same order of magnitude 

as that at v = ^, This also means that the effective diameter 

of the molecules change very little in the shock region for 
•1 

v> ^ and the molecules hehave very nearly as rigid spheres for 
strong shocks. The work on B-G--K model by In ant ha, say an am and 
JTarasimha^ also suggests the choice of Xg. 

Quite often a mean free path based on a collisional pro- 
cess in the region of interest has been used, c.g,. in the 
2 ,tom-atoni collisional ionization dominated region, the atoin- 
itom collisional ionization nean free path^"^ and so on. Chuhh"*^ 
.n his studies has used the atom-electron inelastic collision 
lean free path for nondimensionalising the relaxation zone, 
his length tends to infinity in the region where the gas is 
ully ionized since the atom density is zero. Also, for certain 


1 
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iq 

gases ©.g. ^enon ' , tihis lengtiTi "becones 1;oo large because of 
the -very small collisional cross-section and hence cannot be 
used as an appropriate scaling length, leverthelesss the choice 
of the length sca,le is not serious beca,use results from one 
length to another con always be transformed without -undue diffi- 
culty. In the present study sti eltogetho-r difierent length scale 
is used.Yfe construct a length scale here ^nich in some sense is 
not a 'mean free path' . We talce the upstream riass flow velo- 
city instead of the upstream therma,! speed of the gas for 
constructing the length. T'his length is defined as follows. 


I 


U. 1/2 

?: (ag-) ■ 


an 


( 4 . 64 ) 


0 aa 

We can show that 1 is proportional to the mean free path 
on the hot side of the gas for large Ma,ch numbers and hence is 
finite for infinite ;'/iach numbers, from reference 35, 


Ag 




9 


iT 


/ r- \ T"! / f-T //-\ \ 


t;1/4 


( 4 . 65 ) 
= 5 , 


/here CT is the effective diaiucter of the molecule, 
'hen. 


f2%n2 A2(5)r(3/2) ""aa 




-2 


( 4 . 66 ) 


ow frori the R-H relations for a normal shock, as 
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Y-«-1 





u 


2 

1 



where y is the racio of specific heats. 


Thus from eq. (4.66), 

X T' V 2 ^ T'/ T' V TT 

2 ^2 h^. 

• • ^ 2 ^ large 

The same length has also been used ■b 3 r Oherai^®. This 
length 1 has all the advantages of the dovmstreani mean free 
path hut still denotes the upstream flow which is ]mown before 
hand both in experimental and theoretical investigations. 

4 . 6 Ho ndLncnsi on a,! Eq uation s ; 

The equations described in the Section (4.4) have been 
nondimensionalized using the parameters discussed in the Sec- 
tion (4.5) and the resulting equations are given below. 


4.6,1 Atom Equations; 


d 


r- 

dx q 


■4- ^ 


1 ^aa 

TS’ 17'-^ 
ae 


. 3/2 =„" 5/2 

. (1-D9 R 


cl 


]'n1 . N 

(-2^ / -i- ) 

^P1 ^ ^‘a 


J 5 - q h 


+ (l-Dj \ ) 


‘1 ■*' "a 


. If, F?/2 e, (is., 

^ ^ fPje 


(4.67) 
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where 


and 


D. 


I). 


V2 fae, «3/2 


i2' 




m -m 


iri 


' Pi 


a.a 


4„y2 y 




(4-. 68) 


[ Z H. U.(u2 + I p ) ] 


dx 


D 3 ' D 


U^i E 
«] 


D 

1 ..2 ^ae V2 

1 AgO) Z ’4 . f2 (g.^-) 


N.p 

e D^oe 


P-i B^, 1/2 

t;; ’ W’ ^ ^e;;) 


'^e ^36 ' f-oe 

Pn *^”11 i 

p J.'j ^ P ) «1 J, 

' Pdi ^^P3i 2 ^ 


5/2 


„„.iL 4 - J. AT AT o' 

■^9 ^ifF s Q j ITr s \rfi ' ^ 2 jf^je 


h *^1. 

■'2 TTjJ’ 

’2 f/’ fp/e"’ 'VTjF 

•(tr^ +1 ep(2 +£g)]] 




5/2 


Dp Q 

■/. h 

3 1 e 


(4.69) 


4.6,2 Electron Eq.u at_in ns : 


d 


p Pe't] = - -F ft 

dx dx 3 'J 


(4.70) 


~ [HeCEe ^ T Pe)]=h 
TP,- 


-E ITg •^ Ej Hg Z HjUjg+B^A^(2) 

D P 


1 e 


le 


[®^P(-trL/Pel)- ^ 


Tc ’ ^ie -0^4:. / '*^ei 


U 


ei 


Erf 


(4.71) 
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where 


f 2 (. 

'■"ae s 

1/2 

A^( 5 ) 


aa 



2 J .2 

(?ie 

2 

■) 

ra 1 / 

( a) 

Vtx 



■^aa 

Pe^'e 

C\J CD 

-I 

Pe)] = 


1 


(4.72) 


J 


"D f -I C‘ ^ ■ T 

* ^2 Mrs*''"”’ "0™“’ vri”’ *’ ’ ' 

VPej ^ej VPeo 


r% h, ^(2)17^1^ 


p Pe ^\e 

• -‘-7, \\Ta" -i 15"--3 sro” 


/ e ' « xe : r \ ! -n-- 

3 M/Tie’ Pie’ ^ ^ 1 


* Pe 


exp (^4,) "A) 


^^ 3/2 U. g 

■ ^ “3 Pje wf-) 


'36 'ih'je 


(4.73) 


whe re , 


4.6.3 Pl_a_^^ Equa t ions ; 

-Sr [I + (1-e) iTA + en;tfj = o, 

dx 3 'J 

v/hich can he written using eq. (4.59) as, 


- 4 - [E SXj + H^uj = 0 

dx D 


(4.74) 
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Eqn. (4.56) after substituting for (n -n.) from (4.62) and 

0 1 

using eq. (4.63)s reduces to, 


dx Lj 3 '' 3 23 

a .2 


r, ae. ' 

4 “0 ”'2" 


= 0 


(4.75) 


and eq, (4.61) reduces to 


d 

dS 


S + !■ P^) + (1-S) + I p, ) 


e i^eUp.(u^ + I Pp + Si) 


= 0 


(4.76) 


•e^e^'^e " 2 Pe 
4.7 _Bqundag 7 _ ..Qo^d.j.t_io_ns ; 

TTne boundary conditions go'veming tbe various parame- 


ters in tbe 

distribution 

function 

z - <=° ; 

1 , 

>5 - 0. 

and 



X ->- + 005* 

N-j 0 5 

S 2 0 9 


Ue , Ui r 

(-U 3 ); 1 


(4.77b) 


4 . 8 Pla sma Equ ations : 

Equations (4.74) to (4.76) can be integrated and tbe 
boundary conditions (4.77a) if introduced, yields, 


Z 

3 


“Pd ^ 




U- 


(4.78) 
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s ^ (1..) + 1 p ) ^ 


'sC'ff + i Pe) 


2£ f 2 ""aa i^*2 

4 " --0 ^ 


o ^ 

:u; I- 4; 


(4.79) 


f 7’7(’^5 I PP " V4 ! (i-s)(t? + I- p,.) 


S'l) f 


1 2' I'l' 


■J'lCtff + 1) 


(4.80) 


Equations (4.78) to (4.80) represent the plasina equations of 
conservation of mass, momentum and energy within the shockwave, 
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iVaSTiiOD OE SOLUTION MID PHTSICAL PARMETERS 

5.1 G^^eral; 

In this chapter the nondimensional equations derived 
in Chapter 4 will he simplified to describe the flow in the 
downstream equilibrium region by introducing appropriate 
bound8,ry conditions and the resulting equ,ations will be 
numerically solved for the junp quantities across the shock, 
I'he structure equations also will be solved to get the shock 
structure profiles for a given set of initial conditions till 
the final equilibriun values are reached. Various physical 
parameters of interest are computed. 

5.2 E q uation _s _ Bepr e sent ing the J ump Pond it ions s 

Introducing the boundary condition (4.77b) in the 
equations (4.78) to (4.80), denoting the values of 
R-, "P '7 for X + =°5 by h*, H*, H*, "f''' and using, 

Pe = P*/® ’ Pi " FV(i-e) 

yields, 

(R* + IT*) U* = (5.1) 

-I- ^ p*) 4- (1-e) R*(IJ*^ + ^) 

+ S lTg(U*'^ + p ) ~ (U^ + •^) 


(5.2) 
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:^2 ^ 

1 B *' 



£ 

( 11*2 

c 4 - 

D "* 

CO 

CD 






-r £ (Tj-- + .^t’ ?" + gj)j= *^1)“ (5.3) 

Introducing a, the degree of ionization, defined hy 

a = = ^ ~ J-'^eVCh'^+H'^) (5. 

equations (5.1) to (5.3) can he rewritten as 

Ij* U* = tJF(1-a'') (5.5) 

11 + %- (1+h)l = (1-a )(U? + wj (5.6) 


1 /(M +t ) , 
e' ^ e a-' ^ 


</(he+ir) (5.4) 


(1+tf)] = (1-a'')(u2 + -1; 


l4''~ [11*2 ^ 5 ^ = (l~a*)‘tj^(tJ^ + -1) 


* „2 


f2 5- 

H ^ 2- 


I'he Saha equation for equilibriun, is 

,..yi = 1. a! p*V2 exp (-..n-) 

1-a-' ^o ^1 kly:* 


. 5 , 1 ) 


(5.8) 


where , 


2g, 2nra kiy 3/2 

■"T / 0 1 y 

. V — — g— ^ 
h 


■Jiliminating 'p* between eqtiations (3.6) and (5.7) a^id substi- 
tuting from eq. (5.5), yields. 


I RS , o- -4)^ 


■f5’ ^2' 


where is defined by. 


5 2 

’F” '^1 
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The positive sign behind the radical sign in eq. (5.9) 
has been oraitted as it corresponds to a, trivial solution when 
= 0. From eqs, (5.5) and (5.9), it ccn be written that. 


f- (l~a''')/r(|- 


■I ^ /.'‘-i ^ ^2 32 a'' . I s 

’‘ 2 )-'* O'- ‘* 2 ^ "i'5' ■'■'2’ 

■-I M. '2 


(5.10) 

From equations ( 5 . 5)5 (5.6) and (5.7), 'p* can be solved. 
Using eq. ( 5 . 9)5 "be written as 


P'' = 


2 (1-fa*) 

- /( 1 - 


1 32 a* 

e- I ^ 7v/r f- 

"‘1 


4 - 

V~“ 

w2) : 
"1^ i 

\}'3 

,.1 .7 

Mgj - 

I 

klT,' 

* a 


(5.11) 


Equations ( 5 . 9)5 (5.10) and (5.11) give e^rplicit expressions 
for the velocity, density and temperature across the shock 
in terms of the initial l[ach number and equilibrium degree 
of ionization. 

Equation (5.8) can be written using eq. (5.10), as 

- *2 ^2 rr Jr-, r~ ir -5 


= 5 n e*V2 55 , .J. ^ 

8 n„ P Ih 2> 


(1- 1 )P + f ( I )1 exp (-lAiy*) 

M-1 }\x^ 1 — ^ 

' ^ ( 5 . 12 ) 


Equations (5.9) to (5.12) constitute the expressions for the 
jump conditions across an ionizing shockwave. 
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JiJquatioiis (5.11) and (5»12) have heeii simultaneously 
solved for p''"' and a" numerically usin 3 , I'evreon-Raphson method 
with double precision on I'Bil 7044 Comrjuter for argon , krypton 
and xenon for initial pressures of 0.1 to 100 torrs with 
initial temperature = 300°Iv and Mach numbers 'ranging from 
10 to 40. Ihe velocity ratio and the density ratio have been 
calculated from eqs. (5.9) and (5.10). Ihe Mach number for 
complete ionization of atoms has also been computec’ . The 
various physical parameters that will be discussed in Section 
(5.4) have been computed and the results presented in tables 
(5.1 - 5.6). 

5 . 3 Shock _St r'ucture_ Pro f ; 

The set of equations (4.67) to ( 4 . 76 ) and (4,62) can be 
solved simultaneously subject to the initial conditions (4.77a.). 
Equations (4.62) ?;ill be dropped and the qua,si~ueutrality assump- 
tion iT = N. is made eveiywheie. Using equation (4.60), this 
yields IJ„ = Ut for all x. Since m.'.^ m and the atom-ion 

vi? _L JL Cv 

collision cross section is very high, the atom and ion tempera- 
tures can be assumed to be same aiid hence everywhere T =1)”.. 

Also it can be ass'umed that Ll = This gives, 

and 




T. 

1 


(5,13) 
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With these approximations, the moment-um conservation 

for plasma is aiaproximately satisfied and hence equation (4,75) 

can he ignored, further li << m , << and hence the 

G a, 6 ‘^6 

kinetic energy of mass motion of electrons can he neglected 
compared to its thermal energy. With these assumptions, we 
have the follov/inf equations. 


h = 1 B. V 


dx 


4 

<) 




'B.e 


where 


and 


R. 


R. 


Jif 




■ e-^^Pje ^1Tp^e’ fp^e 
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3/2 ^.- 3/2 


2 “ Pa 




..- 3/2 .e. 


( Wl ) 


Dg Pe 

... 1 

2heUi dx 


e^ H 


a 


(BpP 


e^e 


(5.14) 


( 5 . 15 ) 


( 5 . 16 ) 


and 


5,.. .A. 

dx 


Ped'Q 


-2'lU 


e e' 


+ BjHp £ il^ 

J 



H-6pg)+2tJ,(A-U, 


D ' D e 


t 21^142) (Pi-epg) 


le 


1 -2 , S .-^2 B. g„g- 

l “181 f ®j he 
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■je *^^je. 
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which alter using equation (5.16) can he rev^ritten as 


5 


'■ U 
e e 


dx 




U.ID3 Ng z n. 

3 

I'T.l'I 







+ ^D^A^(2) (R^-£Pg) 


'le 


, 3/2 


U.. 


- ffe - -I- We i- (Ve) (5.17) 


In eq. (5.17) the last tern on the right hand side 
represents the ai/erage energy the electron gas has to spare 
for e-very newly bom electron to raise it to the local electron 
gas temperature. If the newly bom electrons are asstiined to 
have already attained the local temperature of the electron gas, 
the last term drops out. Hoffert and. Lin-^ hawe shorn that the 
last term in eq, (5.17) can be ignored in -view of the argument 
that this terii getf-j, cancelled with a similar term arising from 
the time derivative term of Boltmann equation on left hand 
side. In the present work this tern has been omitted. 
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The - moment equation for axioms is. 


xi- & i Pi> ’^i ? 


dx j 




+ ^fz ( a£) 
*^aa 

aa 


B»», 1/2 1/2 


U, 


P 


6. U. 

r ^ -n i 


Vjpje ^1%^’ jfz 


u, p. u.. . 

'l‘'3 ^31 ■?-/ 


dJ 1 .^aa . -^5/2 j, 4°^) 

V D khk Pjk ^^dlc 

u. p. U,, 


ae k 
,•5/2 


J- ‘ TO* }\ o^/ T? ^.1. 

? "e*'3 Pje ^aXrPje’ Pje’ iTPje ’ 


•3' iiple (d'i •<• i Ph(2+ ~°-?.) 


,2 


P. 


(5.18) 


The plasma conservation eouationscan he written from 
eqLuations (4.74) and (4.76) after introducing the boundary 
conditions (4.77a)5 as 


and 


2 N .U. = 
D 


Z i'i4jj(U^ + I P;5) = ‘■CV'Jl + 1) 


(5.19) 


leU, 


Je Pe +g 1 ) ( 5 . 20 ) 

The set of equations (5.14) to (5.20) have been solved 
numerically for the variables 11^, p^, k, , I 2 and using 
Runge-Kutta-Oill’ s method on the IBM 7044 Computer. To initiate 
the integration process a value of 0.99 for and the corres- 
ponding value of I 2 given by the inert gas shock profiles from 
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Mott-iSmith' s solution hauc iDeen used and this initial point has 
heen aroitrarily taken as the origin lor the shock coordinate. 

It is to be recalled that the position of a normal shockwave 
in one dimensional flovj is indeterminate and hence the choice 
of the origin ior the shock coordinate is immaterial, Ihe struc- 
ture profile for negative x will be little different from that 
of inert gas profiles as the degree of ionization is negligibly 
small in this region if the effect of precursor is neglected. 

In the present work the precursor effect is neglected and an 

6 

initial value of 10“ for the degree of ionization at the start- 
ing point has been used. This value ha.s little influence .on 
the structure profiles as reported by Chubb and also checked 
in the present studies for the range of ^ ach numbers studied. 

The initial electron temperature has been taken to be the same 
as the initial atom temperature in the upstream enuilibrium 
region. The electron temperature however-j increases very fast 
initially and remains constant over a large part of the shock,, 
and the value of initial electron temperature assumed affects 
the shock structure negligibly^ . However, the position of the 
onset point varies considerably. 

The following definitions of density, velocity and tem- 
peratures are used. 

n = /f^, U = ™/ fc.^^, |■kT=^m/fc^^, (5.21) 

The parameters of interest namely, temperatures of atom, 
ion and electrons, their velocity, pressure, densities, degree of 
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ionization, electrical conduct ivity^ the electric field and 
potential and refractive index have been conouted at every 
point within the shockwave, ihe integration is continued 
till 11^5 reach values less than 10“"^’. At this stage the 
various quantities v±ll have reached their equilibrium values 
to the third digit accuracy. These results are discussed in 
Chapter 6. 

5*4 Phy si cal Par am e ter 3 : 

1 ) Pressure: 

The pressure in an ionized gas can be ?</ritten as 

p = 2 p. + Ap 

i 

where p. is the partial pressure due to i-th species and Ap 
J- “ c 

is the pressure correction tom arising from the coulomb 

interaction between charged particles. Ap^ is negligible 

2 

since the mean coulomb potential energj^ e /n^ 9 is much less 

3 1/3 

than the mean thermal energy 4 kl , where r = n“ ■ , within 

w tr w O 

the range of present intere&l. Substituting 

p . = n .kT . , 

■*=^1 11 ’ 

the ratio of the local pressure to the free-stream pressure 
can be written as, 
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where p is the dinensionless mass density. At equilibriim, 
Tg - and hence eq. (5.22) reduces to 

P = p (1 + a*)- (5.23) 

2 ) Do TO s t_r eajn . M ac h_ Ilumh ex ; 


The Blach number in the doTOstream equilibrium region 


is 


Mg 


P'4. 

u 


a 


where 'a' is the speed of sound in the equilibriuin region. An 
expression for 'a' can be obtained from the isentropic equa- 
tion, the equation of state and the internal energy expression 
for an ionized gas. The equation of state and the internal 
energy are given by 

p = f^ (p, T) = p 1 (1 + a*(p, T)) (5.24) 

B, 

and 

e = fg (p, T) = I T (1 + a*(p, T)) + -5- a*(p, T) 

'"a ‘ 3. 

(5.25) 

where ^ is the internal energy per unit mass of the gas. 


Eqs. (5.24), (5.25) and the isentropic expressions in 


the differential forms are; 


dp = 


de 


af^ 

~ 

afg 

“BF 


d p + 


dp + 


a f ^ 
TT 
ef 


dT 


FT' 


- dT 


(5.26) 

(5.27) 


d& + pd (J-) = TdS s 0 

-p 

From eqs. (5,26) - (5.28) it can be written that. 


(5.28) 
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T5-’ + ( 


0 f 2 £ f . /S I' 


op 


^■) 


L jL- 




(5.29) 


*> (5.30) 


Using expressions for f-] and fg from eqs. (5,24) and (5.25) 
and for a from eq, (5.8), it can be shown that 

3(4+5a*-5a^'2)e2H.4o,*(-i^Q^*)(3Q+^^ 

where , 

9 = (l/kl')“^ 

"We shall define an adiabatic indez Y for an ionized gas 
by the expression, 

2 


Y 


a' 


P/P 


a 


m 

a 


d+a^'*-) 


which reduces to 

Y = <( 3 ( 4+9a* -5a*^ ) 9^ _± 4.a*( ) ( 59+1 ) 1 

[5(4+5a'"-5a*2)e2+4a*( l_ad (39+1 ) J ( 1+a*) 

Y 5/3 s as rd - 0, 1. 

5 • El e ctrical _ Condu ct ivi ty ; 


The contribution to electricsil conductivity . of an 

ionized gas arises from the close encounters between the 

electrons and neutral particles, and from the long range 

44 

interactions between electrons and ions. Lin et al. have 
pointed out that for a high degree of ionization (a > 10 
for Ar) , the diffusion of electrons is primarily limited by 
long range coulomb interaction of the gas and the contribution 
arising from the electron interaction with neutral particles 
is negligible. 
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Often the Chapmann-Enskog method of solving the Maxwell- 
Bolt2rmaxii'i equation has heen used to get tne electrical conduc- 
tivity of completely ionized gases to various degrees of approxi. 
mat ion. Spitzer et al. ' solved the problem in a different way 
by considering the effect of electron ion encounters as a prob- 
lem of diffusion in the velocity space and their results agree 
closely with that of Chapraann and Cio-wling. 


7. 1 

According to Spitzer and Harm 
tivity of a fully ionized gas is 


CT = 


0.591 (kTg) 


3A 


ill 


172' -2' 


e In ( 


the electrical conduc- 


(5.32) 


where 5 



klg 1/2 

= ( ^ the iJebye shielding length, and 

Sang e 

the landau distance. 


Equation (5.32) has a relative inancuracy of the order 
of In arising from the uncertainty in the choice of the 

cut-off distance. 


In the present case a > 10“"" over a large part of the 
shock and the electrical conauctivitj? calculated from equation 
(5.32) is fairly accurate. 

Re __^_active _Indexs 

For a partially ionized gas composed of ground state 
atoms, ions and free electrons the expression for the refractive 
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index can te written 

p = 1^-2^ (5^n^ +E.n.) - ^ ( |a)" (5.53) 

where w = the plasma frequency, e /m ) j w = impressed 

ir C 

frequency; g , g. are the polarizabilities of the groimd state 

cl 1 

atom and ion respectively, and y the refractive index. The 
expression (5.33) is valid when the plasma, frequency is much 
smaller than the applied frequency and v/hen the number density 
of the excited atoms is small compared to the ground state atom 
mmiber density. The refractive index has been calculated for 
Sodium yellow line ( ^ = 5890 i} . 


5 . E le ct r i c P o t e nt i al ; 

The electric potential due to the field E can be written as 

X 

Y = - f E dx 


o 

which can be stiitably nondiroensionalised and vnritten as 

'KT 



e / E dx 
0 


V 


(5.34). 



CHIPTER 6 


Discussions AsW COUCIUSIOITS 


6.1 Eq,uilil)ri'can Properties: 

The dovmstreaiii equilihriiiia properties computed 
ntruerically for Mach numhers from 10 to 40 for argon, krypton 
and xenon for initial pressures of 0.1, 1,0^ 10 and 100 forrs 
are presented in table s 5.1 to 5.6 for an initial temperature 
of 500°K. Results for argon are also plotted in Pigs. (6. 16, 1) 
to (6,21.1) for different initial pressures and temperatures. 
It is seen that the equilibrium degree of ionization is higher 
for higher initial Mach numbers and temperatures and lower 
initial pressures. The equilibrium temperature ratio increases 
with Mach number but takes much lower "values than in the 
inert gas case. For low initial pressures due to the high 
degree of equilibrium ionization the downstream temperature 
will be lov/er than thaf for higher pressures. The density 
ws curves exhibit a peak where the density attains a 
maximum which can be explained by eq. (5.9) where U /u^ 
is the reciprocal of the density ratio. The second tcna 
under the radical sign in eq. (5,9) causes this type of 
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behaviour, A't low Mach numbers a increases faster than 
so that increases whereby an increase in the Mach 

number results in an increase of the density ratio. At higher 
Mach numbers a* v/ill be close to its upper bound value of 1 
and increases very slowly so that a*/M^ decreases with increase 
of causing the density ratio to decrease, Por very large 
Mach numbers when a*/M^-*0 the density ratio reaches a value 
of 4 corresponding to that for a monatomic gas as may be expected. 
The Mach number for which the equilibrium density is maximum 
increases with increase of upstream pressure. The equivalent 
Y <1,67 (its value for inert monatomic gas) over a range of 
Mach numbers depending upon the initial pressure. At low Mach 
numbers when ionization is insignificant its value is y = 1.67? 
again for high Mach numbers when the gas is fully ionized and 
can be considered as a mixture of monatomic gases of different 
species it attains the value y = 1,67. For low initial 
pressures the variation of y is veiy sudden exhibiting a bucket 
type of behaviour (Fig. 6.20.1) and as the initial pressure 

increases the variations become less steep. The lowest value 


** In the following the words ’high’ and ’low* are used in a 
restricted sense. low corresponds to the region towards 
the lower end of the range of parameter studied and high 
corresponds to the region towards the upper end of it. 
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of "1 reached depends upoii dhe initial pressure and increases 
with increase of the initial pressure. In fig. (6. 21 . 1 ) it is 
seen that the curT/e of the Mach number for > 93 % versus 

upstream pressure for different initial temperatures shows that 
the Mach number increases with pressures the variation becoming 
less as the initial temperature increases. Furthers the Mach 
number for > 99 % increases with decrease of initial 

temperature because as may be expected an increase in either 
the initial temperature or the Initial Mach number is required 
to produce the high temperature behind the shockwave. 

6.2 Shock Structure Profiles; 

The structure profiles for various parameters of 
interest are plotted in Figs. (6. 1.1) to (6.15.1) for argon, 
kiypton and xenon for different initial conditions. On the 
whole the structure profiles show a similar behaviour for 
different initial conditions. Over almost the entire part of 
the gasdynamic shock region the structure profiles are 
independent of the initial pressure and tempera,ture and are 
also independent of the gas. The length of the relaxation 
zone, however, does depend upon the initial conditions like 
pressure, temperature and also the gas. This is pronounced at 
high initial Mach numbers. The intennediate zones do not depend 
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upon "the initial conditions for low Mach numhers whereas for 
higher Mach numbers they do show a weak dependence. 

Curwes in Figs. (6.6.1) to (6.6.3) show typical 

variations of the weight functions N-] » ^^2 ^3 ^sed in the 

trimodal distribution function for the atom gas. In the region 
of the gasdynamic shock, decreases monotonically to zero and 
increases monotonically from zero as in the case of Mott- 
Sniith's study. Towards the end of the gasdynamic shock region 
IT 2 reaches a nazimum and starts decreasing when will have 

started building up. reaches zero g.uite early, particularly, 
for low Mach numberst;. ITg and IT^ will be active throughout the 
relaxation zone. Towards the end of the relaxation zone IT 2 0 
while IT^ approaches the equilibrium atom density value, the 
downstream equilibrium region being characterized by alone. 
"When either alone or IT^ alone is present alongwith IT 2 the 
gasdynamic shock region and the relaxation zone are clearly 
separable, at least, conceptually. For high Mach numbers, 
particularly with the choice of a larger value of (as used 

by earlier investigators) such a separation does not seem to 
be a good approximation and a trimodal representation of the 
present type is perhaps better. 
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The atom and ion tempe rat urea in or ease very rapidly 
within the gasdynamic shock and starts drooping in the 
relaxation zone with large gradients towards the end. The 
maximnm atom and ion temperatures reached will he less than 
the inert gas equilibrium value. The electron temperature 
rises rapidly from its upstream equilibrium value to around 

its downstream equilibrium value remaining almost constant 
over the rest of the flow region where the electron gas loses 
energy in inelastic collisions at just the same rate as it 
gains energy through elastic collisions with heavy particles. 
The initial rise in the electron temperature results from a 
continuous supply of energy through elastic collisions from 
the heavy particles in the absence of inelastic collisions. 

This energy transfer rate depends upon the mass ratio 6 and 
the difference between atom and electron temperatures. For 
high Mach numbers the heavy particle temperatures attained 
within the gasdynamic shock are very high- Consequently the 
electron temperature rises rapidly and reaches a constant 
value almost by the end of the gasdynamic shock region. But 
for low Mach numbers, because of the relatively lower heavy 
particle temperatures developed within the shock the region 
of electron temperature rise extends much beyond the atom-atom 
shock as can he seen from Big. (6.2. 1) . Bor argon due to its 
relatively larger mass ratio, the electron temperature rises 
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fas't and at^fcains an almost; cons'ban't value, whereas for xenon 
due "fco ifs much lower mass radio dhe elecdron demperadure 
nod only rises slowly hud also addains a lower consdand value. 
Over a large part of dhe relaxadion zone dhe elecdron 

dempera.dure rises veiy slowly and reaches dhe adorn demperadure 
dowards dhe end of the relaxation zone and finally hodh dhe 
demperadures dogedher go do dhe eciuilihrium demperadure. 
Similar qualidadive behaviour of dhe elecdron demperadure has 
been reported by earlier invesdigadions also neglecding 
radiadion cooling effecds. for M = 40, nearly half dhe 
elecdron demperadure rise dakes place in a narrow region 
dowards dhe end of dhe relaxadion zone. This is possibly 
due do dhe very low inelasdic energy loss in dhis region 
since recombinadion almosd balances dhe ionizadion rade. Jor 
lower Mach numbers however dhe recombinadion rades were 

found do be less dhen dhe ionizadion rade when dhe dwo 

demperadures ged equalised. For low Mach numbers id was 
nn 

found dhad recombinadion plays a minor role over dhe endire 
relaxadion zone and, hence, could be neglecded. For high 
Mach numbers, however, dhe recombinadion is impordand in 
dedermining dhe elecdron demperadure profiles and, perhaps, 
radiadion cooling also may have considerable influence . 

This aspecd however needs furdher experimendal invesdigadion. 
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The density profiles are practically unaffected 
hy the presence of ionization in the gasdsmamic shock hut 
for a slight rise in density towards the end of the gasdynamic 
shock noticeable only for high Mach numbers. Apart from 
the gasdynamic shock, a large part of the density changes 

occur only towards the end of the relaxation zone where a 
major part of the eq.uilibri'um electron density is built up. 

In this region both the total density and ionization degree 
increase to their equilibrium value rapidly. 

The computed electrical conductivity a has been 
plotted for only that part of the shock where a > 10” since, 
as already mentioned, the present calculations of the 
electrical conductivity are reliable only for a degree of 
ionization greater than 10"^. Qualitatively the behaviour 
of ^ is mainly decided by the shape of the electron 
temperature profiles; the electron density influences 
through a logarithmic term. In the earli,er parts of the 
shock cr increases very fast due to a similar increase in the 
electron temperature. In a large part of the relaxation zone 
where the electron temperature changes veiy little the 
electrical conductivity continues to increase slowly due 
to an increase in the electron density. Towards the end 
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of the relaxation zone once again a fall in electron 
temperature, if any, influences the electrical conductivity 
so that the o profiles show a peak towards the end before 
attaining the equilibrium value, from Pigs. (6.2.4) and 

(6.4.4) it can be observed that for initial pressures of 
1.0 and 10 torrs. the fall in electron temperature is 
sufficient enough to produce a peak in the o cxirves. But 
for an initial pressia-e of 0.1 torr. since the electron 
temperature gradient is small no peak is seen in the a curves. 
Por M = 40, asso ciated with a steep rise in electron 
temperature towards the end of the relaxation zone there is 
also seen a corresponding rapid increase in the electrical 
conduct ivity. 

The refractive index(u) profiles within the shock 

37 

for argon have been obtained using a value of 0.67 

In the earlier region of the shock when the lectron density 
is voiy small the contribution to u arises only from the 
second term of the eq. (5.33) and increases as the density 
increases. Towards the downstream end where the electron 
density is sufficient enough so that the third term of 
eq, (5.33) becomes comparable to the second, p ceases to 
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increase and starts decreasing furtlaer downstream. The 

refractive index profiles are in qualitative agreement with 

those obtained from the observed fringe shifts in the 
37 77 

experiments. ’ 

^ongst the other parameters of interest the profiles 
of electric field and the electric potential are also presented. 
Thus, fig. (6, 1 1 , 1 )shows a typical variation of electric field 
with two peaks. At first, because of the large electron 

temperature gradients present a themal diffusion of electrons 
takes place giving rise to a steep rise in the electric field. 

As the electron temperature gradient decreases to a very low 
value the electric field also decreases and then slowly 
increases because of the electron diffusion resulting from 
the electron density gradient, over a large part of the 
relaxa,tiou zone. Towards the end of the relaxation zone 
because of veiy large electron density gradients causing 
strong electron diffusion the electric field again increases 
rapidly and ultimately decreases to zero at the downstream end 
where no gradients are present. The electric potential 

(figs, 6.7.1 to 6.7.3) is positive every where which indicates 
that the electrons move upstream relative to the ions to 
smoothout the gradients. The other parameters of interest 
namely the pressure, velocity etc. are also presented and 
need no explanation. 
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Ctase't-ciegree of Ionizati on; 

The onset-degree of ionization is of interest as it 
helps to identify the region of atoE-electron collisionaa 
ionization dominance. The region downstream of the onset 
point is practically independent of the priming mechanisms. 

The onset-a and the onset point for different shock Mach 
numbers are given in tables (6,1) to (6.6). The onset-a is 
also plotted against Mach number in Pig. (6,22.1). It can be 
seen that the onset degree of ionization increases with increase 
of Mach number and seems to reach an asymptotic value for 
large Mach numbers. However, at such high Mach numbers it 
is doubtful if the atom-atom collisional ionization is the 
appropriate initiating mechanism. The onset-a strongly 
depends on the value of has a higher value for kiypton 

and a relatively lower value for xenon. The onset degree of 
ionizolion has a slight dependence on the initial pressure. 

6,4- Shock Thickness: 

The shock thickness and the thickness parameter 
have been computed for each case and are presented in tables 
6.1 to 6,6 for the three gases. Since collisional ionization 
is predominantly binary the 2 ?elaxation time will exhibit an 
inverse initial pressure dependence. Conseq.uently , t is 
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independent of the initial pressure and has been recommended as 
a suitable parameter for comparison purposes. The shock 
thickness 6 decreases with increase of Mach number. Since the 
nondimens tonal shock thickness increases i7ery slowly with 
pressure, the actual shock thickness decreases almost inversely 
with pressure. The shock thickness also depends upon the 
excitation and ionization potentials and the cross section 
for the atom-atom collisional ionization. It increases with 
decrease of the atom-atom collisional ionization cross section. 

The shock thickness is less sensitive to because the rate 
of atom-electron collision ionization is decided by the rate 
of energy transfer to the electron gas from ion gas, which in 
turn depends upon the mass ratio 6. It may be seen that krypton 
has the smallest shock thickness and xenon has the largest 
amongst the three gases studied. The smaller shock thickness 
in krypton is due to its lower excitation and ionization 
potentials and also its larger atom-atom collisional 
excitation cross section. Since argon has a slightly lower 
walue of this cross section and comparatively larger 
activation energies it has shock thicknesses larger than those 
of krypton gas. On the other hand though xenon has the 
lowest values of activation potentials compared to argon and 
krypton, its atom-atom collisional excitation cross section 
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is nearly 10 times smaller. This in addition to its relatively 
larger atomic weight causes much broader shocks in xenon 
compared to the other two gases. 

The thickness parameter t is plotted against 

Mach number for the three gases in Pig. (6.23.1). The results 

of Chubb and also the results from the present investigation 

for argon with 1.2 x 10 cm /eV have been included for 

comparison. The qualitative behaviour is in general agreement 

with Chubb’s results and the quantitative discrepancy is 

mainly due to a different value of chosen. The curve for 

argon with C^^= 1.2 x 10“"^^ cm^/eY lies very close to the 

corresponding curve of Chubb but has a faster asymptotic 

variation with Mach number than the other. In Pig. (6.24.1) 

the thickness parameter p^ x is plotted against the reciprocal 

of the maximum atom temperature within the shock to compare 

the present results with earlier theories and experiments. 

77 

Some of the results of Wong and Bershader , and Petschek and 

Byron^^ are also reproduced here. The results of Hoffert and 

36 5Q 

lien and also Morgan and Morrison have been given. Results 

of the present investigation seem to show a close agreement 

with the measurements of Wong and Bershader at lower Mach 

numbers. Lack of availability of cxperiment.-xL result* in the 

high Mach number range precludes any better comparison than 
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this. The high Mach nimher results of Hoffert and lien, and 
those of Petschek and Byron for higher degree of purity 
extrapolated to larger Mach nunhers, seem to agree with the 
present studies. 

6,5 Gon elusio ns ; 

Por large Mach numbers (M>-^30) it is necessary to 
consider the gasdynamic shock and the relaxation zone 
simultaneously since it is difficult to assign a definite 
plane of demarcation between the two regions. It is still 
worse when a relatively larger value of the atom-atom 
collisional excitation cross section is used for which there 
does not exist an intermediate plateau in the shock structure 
profiles. A trimodal distribution function model of the atom 
gas on the other hand overcomes this difficulty and the entire 
region of nonequilibrium can be studied uniformly without 
resorting to identifying the two regions, Howeverj no 
justification for a trimodal distribution postulate has been 
given. 

The good agreement of the results of the present 
work with those of Wong and Bershader is mainly attributable 
to the better choice of the atom-atom collisional excitation 
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cross section used. The shock thickness and the onset point 
depend strongly upon the -value of this cross section -whereas 

they are less sensiti-ve to the choice of the aton-electron 
collisional excitation cross section, itirther a steady state 
electron energy equation seems to he a good approximation for 
a large part of the relaxation zone and the electron temperature 
variations are significant only in the gasdynamic shock and 
to-wards the end of the relaxation zone. 

The effect of radiation on initial ionization and 
cooling of the downstream eq-uilihri-um region a,long with 
absorption and emission of the radiation hy the gas is to he 
included for a more meaningful study. More reliable experimental 
cross sections for the photoprocesses still remain for 
investigation before a general study of this type can be 
carried through. The effect of second ionization is to be 
included in large Mach n-umber studies. 
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Table 6,1 

Gbaract eristics of Ionizing Sliock 




Argon : p 

^ = 0.1 

torr. T 

^ = 300°K 



I ( cms ) 

a~0n set 
3C10^ 

x-Onset 

x10“^ 

6 ( eras) 
x10~^ 

T (secs. 
xIO^' 

) 6/L 

_? 

x10 

15 

0.2753 

0.0919 

1 . 3600 

1.4636 

3.0241 

5.3159 

20 

0.3671 

0.4923 

0.3482 

0.3282 

0.5086 

0.8940 

25 

0.4589 

1.0175 

0.1612 

0.1565 

0.1940 

0.3410 

30 

0.5506 

1.5298 

0.1031 

0.1065 

0.1100 

0.1934 

35 

0.6424 

2.0033 

0.0789 

0.0892 

0.0790 

0.1388 

40 

0.7342 

2.4384 

0.0664 

0.0772 

0.0599 

0. 1052 
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Table 6.2 

Characteristics of Ionizing Shock 




Argon: p_^ = 

1.0 torr. 

Tl = 

300 K 


l( cms) 

a- On set 

X- On set 

6( cms) 

T (secs 

.) 5/L 

o 

I 

xIO 

2 

zlO 

-2 

xIO 

xIO"^ 

xIO^ 

xIO" 

15 

0.2753 

0.0940 

1 . 4000 

1.4831 

3.0645 

5.3869 

20 

0.3671 

0.5173 

0.3630 

0.3449 

0.5342 

0.9396 

25 

0.4589 

1.0790 

0. 1684 

0.1664 

0.2063 

0.3626 

50 

0.5506 

1.6303 

0.1073 

0.1122 

0*1159 

0.2037 

35 

0.6424 

2.1377 

0.0815 

0.0936 

0.0829 

0.1457 

40 

0.7342 

2.6075 

0.0683 

0.0809 

0.0627 

0.1102 
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Table 6,3 

Characteristics of Ionizing Shock 


Argon: = 10,0 torrs, = 300°K 



1( cms) 
2 

z10^ 

a-Onset 

2 

xIO 

z~0nset 

x10”^ 

6 ( cm s) 

X (secs. 
x10^ 

) 5/1 

-2 

x10 

15 

0.2753 

0.0956 

1.4100 

1.5233 

3.1476 

5.533 

20 

0.3671 

0.5469 

0.3805 

0.3682 

0.5706 

1 ,0030 

25 

0.4589 

1.1577 

0.1775 

0. 1796 

0.2226 

0.3913 

30 

0 . 5506 

1.7578 

0.1126 

0. 1211 

0.1251 

0.2199 

35 

0.6424 

2.3111 

0 . 8490 

0. 0998 

0.0883 

0.1553 

40 

0.7342 

2.8083 

0.7060 

0.0859 

0.0666 

0.1170 
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To-ble 6.4 

Glia,rncterif-tics of Ionizing Shock 




Argon; 

aa 

= 1.0 torn. = 300"Ii 

= 1.2x10“^^ cx^/eY 



L( cms) 

a~0nset 

x-Onset 

6 ( cms ) 

T ( secs. 

) 6/L 


x10 

xIO^ 

x10" 

-2 

xIO 

x1o5 

_2 

x10 

15 

0.2753 

0.3283 

1.0620 


0.9017 

1.8631 

3.2750 

20 

0.3670 

1.4235 

0.2355 


0. 1920 

0.2975 

0.5230 

25 

0.4589 

2,7642 

0.1119 


0.0941 

0.1167 

0.2051 

30 

0.5506 

4.0581 

0.0762 


0.0672 

0.0695 

0.1221 

35 

0.6424 

5.2386 

0.0619 


0.0614 

0.0543 

0.0955 

40 

0.7342 

6.3761 

0.0548 


0.0601 

0.0466 

0.0819 



149 


Table 6.5 

Characteristics of lonizin.:? Shock 




Krypton; = 

1 .0 torr. 

T., = 300° 

K 


L( cms) 

a-Onset 

o 

X- On set 

6 ( a ■! s ) 

T ( secs) 

6/1 


xIO 

xio"^ 

x10“^ 

xl0“'' 

xIO^ 

-2 

xIO 

15 

0. 1788 

0.1707 

1.6470 

0.9770 

2.9222 

5.4650 

20 

0.2384 

0.7159 

0.4497 

0.2650 

0.5945 

1.1117 

25 

0.2980 

1.3563 

0.2111 

0.1359 

0.2439 

0.4561 

50 

0.3575 

1.9586 

0.1328 

0.0935 

0.1398 

0.2615 

35 

0.4175 

2.5073 

0.0987 

0.0780 

0.1000 

0.1870 

40 

0.4767 

3.0032 

0.0806 

0.0697 

0.0782 

0.1463 
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Table 6.6 


Characteristics of Ionizing Shock 




Xenon 

■: = 1 

,0 torn. 

T^ = 300°K 


l(cms) 

x10 

a-Onset 

x10^ 

x-Onset 

-2 

x10 

6(cms) 

xIO"^ 

T ( secs. 
x10^ 

) 6/1 

_2 

x10 

15 

0.1377 

0.0521 

2.3240 

1.5118 

5.6634 

10.9820 

20 

0.1835 

0.2092 

0.8192 

0.5080 

1.4274 

2.7680 

25 

0.2980 

0.3982 

0.3340 

0.2854 

0.6414 

0.9577 

30 

0.2753 

0.5775 

0.2662 

0.1975 

0.3700 

0.7175 

55 

0.3212 

0.7429 

0.1904 

0.1558 

0.2501 

0.4850 

40 

0.4767 

0.8918 

0.1222 

0.1382 

0.1942 

0.2900 



APPEKDIS 




All "tlie integrals encountered in tlie evaluation of 
inelastic collision terms come under the general form, 

I = j exp J |;g2.g2;[ g 4 , (^) ^ ^ 

(A-1) 

where all the variables are dimensionless and the domain of inte- 




gration is that part of ? - n spa.ce in which. 


g > Sr 


(A-2) 


where g is defined by, 


g = C - n . 


(J)(C) is the moment function and 'i is the unit vector along 
x-direction. 5 and n have components and ( n^, n^) 

along X, y, z directions. 


Integral (A-1) is evaluated for different 4)(?). 


(a) Tra nsformations: 

Introducing the follov/ing transfomations T^hich simplify 
the terms in the exponent of eq. (A-l), 


and 

„ (u_Oi = S' + (1-P) 


(1-3) 



it can Toe shown that, 




I'" 


Ti-J)""-™- - p(Trp7 


? - n 


Ci - Cl 


and 


(p (0 = p, J, ^). 


Integral (A-1) can he written as 

.2 


/. 


I'',G 


exp 




(A-4) 

(A-5) 

(A-6) 


(g^-gf:) g<i.(U, P, S’, &)J <i& 

(A-7) 


where J is the Jacohian of transformation. Erom eq.s. (A-3) 

1 . 

Introducing the spherical coordinates for ^ and ^ as, 

(E, 0ji, 4j,) 




E 


and 


(A-8) 


^ (G, 0c|5 4-(.) 

the elemental volumes ^ and ^ can he written as 

dF = E^ sin 6j, dE de^, dhp , 

^ sin 0^ dG de^ 

and the relative speed g is given hy 
„2 __ 


(A-9) 


+ q. 2 _ 2^0- cos Bq. 


( 1 - 10 ) 


<N 


since G * i 


G cos 0 


G * 



iii 


'b) Limit_s of Integration: 

From eq. (a~ 10) it car be seen that g is independent of 
and hence bhe constraint (A-2) is a, constraint only on G-. 
Therefore the limits of P are, 

F K--. 0 to 00, ©pv^O to 7X5 (jjp 0 to 2% . (A-11) 

limits of ^ should satisfy (A-2). The minimum value of g is 
g^, so that for every Q^, there is a lov^er limit for & = G (e„) 

0 It 

given by, 

go = + C-2 _ 2 cGo 003 0g . (A-12) 

Cuive represented by eq. (A-12) is shovm in Fig. (A-1) which 
divides the G - Qq domain into two regions, in one g < Sq and 
in the other g > g^. Eq. (A~12) can be solved for to give, 

Gq = ^ cos ©Q + 1/ g^ - sin 0^, . (A-1 3) 

where 5 >_ 0 and for G^ to be positive it can be 

shown that positive sign is the proper choice in (A-1 3). 


It will be shown later that ^ > 0 is not a limitation 
however, and > ? ia satisfied in the range of speeds of 
interest in the present studies. 


Tn Fig. (A-1) is shown the region of integration in G - 9 
plane denoted by R, which can be divided into R^ and R2 as shown, 

so that , r r r 

I = i ^ ( (A— 14) 

R R. R2 


G 



iv 


and the lln.lt a of integration for and are appropriately 
denoted by, 

OO _ 

I = 


R. 


dG 


and 


go+C 


/ 


de(j 


0 




TC 


i = / <1= / 


R2 


de^ 


^o- ^ 


’GO 


where given by, 


S. 


2 


C + G^ - 2 ? G cos 0 


GO 


Finally b q_ varies from 0 to 27 t;. 


(A- 15) 


(A-16) 


° ) Svaluat J-on_ of t he Inte^agJL : 


Introducing eqs. (A-8) and (A-9) into eq. (A-7) and using 
the limits defined earlier, eq. (A-7) can be shown to be equal to. 


OO 

r 

u 2-n; 


I = f 

L f . i 

7 

ii 

o 


0- ^ 

d^l 

f/ F^ 

+ (g^-go) g < 

F^G' 


Rr 


’F 


d0p dG^ d4>p d4ig_ . 


uu 

(A- 17) 


Since F and G are independent vectors, integration over F and 
G can be performed independently. In the following, I is 
evaluated for different moment functions of interest. 



V 


■ I = 1; 

I'rom eq. (ii~. 17) 


r?° 71 

/ [ I 



Stg 

r 

/■ 

-G- 

11 

O 

j 

o 

!i 




r. ]?^ 

L( 


+ C&^-Eq) g sin e^, sin 


dl' dG d9j, deQ_ d',^ dd^ 


(A- 18) 


Inte(iCatin£' over ^ and (|)q_ , and using the results, 

Tt 2% CO 

j' Sin 0 d0 = 2, j dd = 2715 j" exp (-x^) x^ dx = fu/A 

0 0 0 

eq. (A- 18) can be written as 


1 

= 2x5/2 

p5/2 ( 1 . 

f 


exp 

(-g 2 ) 

(g^- 

go) gO^ 




R, 

+R 2 






sin 0 Q. d 0 g_ 

dG . 




(--19) 


of the 

Integral 

Over : 






/ 

exp (-G^ 

X , 2 2 

) (g -So 

) gG'^ sin 

Qg 

de^ 

vT 

dG =- 

00 

f 

exp (~G' 








g 4 

■ C 



% 








X 

.27 

(g^-go) g 

; sin 0 g 

“&1 

dG. 

(A-20) 


0 

Prom eq. (A-IO), 

C G sin 0Q. d0j. = gdg 

and when 0 q, = O; g =?-05 9q^ = x; g = c:+G. 


(A-21) 



VI 


using in eq. (A-20), it can be y/ritten that, 

^ C +& 

/ (s -go) e sin 00 (100 = {g‘-go) g‘ ag, 


0 


which after integration reduces to 
% 


j' g sin 

^ 0 


®e “g = G 


; -G 


sU(252-g2)G2 


+ 1 ,2 1 

+ <= t g^_) 


( 1 - 22 ) 


Introducing eq. (A-22) in eq. (A-20), 


oo 

[ e^q) (-G^)(g^-gQ) gG^ sin 0^ de^ dQ = 2 f exp (-G" 


R, 




G 


G 


+ (2S^-g^) G^ + C 2 ^ g^) 


dG 


which reduces to, 


f 

R. 


exp 


(go+ ?) ' 


^ d * 8go G + g2 + 4) 


+ 2g^ (gQ+2)C + (gQ+2) 


(—23) 


E- valua tion of Integr a l O've r R 2 : 


go+? 


r exp (-G^) (g^-g^) gG^ sin 0^ d0^ dG = J exp (_G^) 

J O' V 


R. 


So-? 


X G 


2 


TC 

j" (g^-g^) 6 sin 0Q_ d0Ql dG (A-24) 
0 GO 


Eroni eqs. (a- 10) and (A-16), 
sin 0Q. d9g. = I 

and when 9^. = Qqo’ g = Sn? ^ + &• 


(-1-25) 


^G 



Using (A-25), it can Ue written that 


2 2 

(u -g") g sin Bp de„ 

O- Ijr 


? + G 


/ 2 2\ 2 
(g-go) g dg 


which after integrating reduces to 


(g^-g^) g sin 0 ^ de^ = y. g5 + ( 2 c 2 „ ^ g2) 


G^ + Q. 


+ (-i _ 4“ g^ IT 4f 

^ 3 °o 15 °o^ 


(A-26) 


Inserting (A-26) in (a- 24) and integrating over G yields, 

/ (-G-^) (g^-gg) gG^ sin Gp dGp dG 

R- e- b- 


% ‘e 


Erf (g^+0 - Erf ^ 4(3_g2)^2 


+ (3-2g^4 + ^ exp|^-(g^_^)2l 


1-4- g + 4- C + 


- n - ^6g^ + (j^ + 40 g2^ ^5 

+ (4g^ + g4 ^ (-^ + 2g2)? + I gj ^_27) 

Adding e^s. (A— 23) and (A— 27)5 and introducing the resulting 
expression in eq. (A~19) yields finally, 

I = 4%3/2 (l.p)V2,^(,, (,_23) 


Where, 



viii 


^0^ ~ Erf (g^+5) _ Erf (g^-C)j 


^ (3-g2)^+ (^ 


f“ “ f Sq) 4‘ I + e 2 cpf-(g 2+52 


X 


(2c +5) cosh 2g c+ (2c2+3)t2-c^ 


and 


+3) sinh 2g C 
C ° - 


■^'1(0 , Sq) - 4(2 + g^) exp (-g^) 

_Ca^_, when _C <_ 0^: 

Tho transformation similar to (..-3) is chosen 


5 (n-29) 

(-^-30) 


as, 


and 


TJ i = E + BCr 


- (U-C)^ = E _ (1_p) 0 


-31) 


so that cq, (^i.— 4) is unaltered hut (*1.— 5) hecoices. 


e = Cl + G 
and hence eq. (.i.~10) is, 

g = + G + 2GC 

= c'^ + - 2G^' 

where , 

C’ (= -C) > 0. 

Therefore eq, („x--7) reduces to 

I - / exp (- 

E,G- P(1-P) 

X (g^-g^) g ^^(U, p, E, J dG ^ (A-53) 

where C ' > 0. 

Thus the result for c < 0 can he obtained hy substituting 
C ’(= -C) in (A-29) and hence (A-29) is valid for ^ ^ 0. 


( 1-32) 



A-II "f" ( 1 ) = ? 


X 


Proi'd eqs. (ii-3) , 

5 U + P cos e _ p G- cos e 


& 


aJid 


i^iCUjPsPsG) - TJ'^+P^ cos 0J, + p^G^ cos^ eQ_ + 2UP cos 9^, 


~ 2UpG cos 9p_ -2PpC:" cos 0^, cos 0g_ . (^^-34) 

Using eq, (2^-34) in (^^-17) and following a similar procedure as 
in it can "be shown that, 


whe re 


. 5/2 , 3/2 


7 t " ( 1 ~ p )^/^ PgCU , p , g ^), 


(A-35) 


P 2 (U,p,C,gQ) = I^Brf (g^+c) - Erf (gQ-4)| 

X j(U^+ p) 4 ^ + 3 Up:^H-[| p^' + (^-g^) (U^ + |)j ? 

+ Up(3-g2) + I (u2 + ^ p) (3-2g2) I „ ^ up (3-2g2) 

+ ^ p (3-2gQ) --3- 


1+ exp (U^+g^) i + '2' 


+ ^(U^ + ^p) + '^p^+-j dp ( 3+4 g^) 

» {■ P^(3+4 g^) -g- j cosh 2g^^ 

.2 


, 2 ^ 1 


+ [(u^ + ^ fj) +'3up C+ I P® (3 + 2go) + I (u^ + 5 
- I up d + I 1 j slnh 2 g„2 . (A- 36) 

Cj, _-J 

Eq. (;i.-36) is valid for 4 ^0, and 

E 2 (U,p, 0 ,g^) = exp (-g^) I p^ gj + (2U^ + | p^ + p) 

+ (4U^ + 2p^ + 2p)J . (A-3T) 



X 


2 . Ela stic ; 


.ii.ll the integrals encotmtered in the e^/aluation of elastic 
collision, terns come under the general fn -m . 


I = 


r 

j 


A ',2 


oxp 


(§- ui) 
“I' 


?^^2 


(a-('J-c)i) 


(IsIlO < 3 ^ 

(-1-38) 


where o depends upon the force law of interaction and the moment 
function, i is the unit vector along x~direction. ^ and x] 

V.-«- 

have components (5^, ^'z^‘ 

a) 1 ra n s f o ma t io n s ^ 


Introducing the transformations (1-3) nwe get, 

(?-U£)2 (ri_(U-c)i)^ p2 2 

(j , (_39) 

g =1,"P,= (-i^-40) 

and 

I’ 

Integral (1-38) can he v/ritten as ! 

I = / exp (- - G-^) °'^(U,p}CjE, ^ (1-42 

E,& i 

w. w^, 

Eurther, spherical coordinates nan he introduced for E and 2^ 
from eqs. (.^..-8) to (1-10). 

h ) limits of Int egr atio n; ' 

Integral (..1-42) is to he evaluated suh;iect to g > 0. It | 

can he le called from Section 1 of Appendix that integration over j 



xi 


^ can "be perfura'ed independent of this constraint. Therefore the 
limits on ^ are the same as given oy' ( ^ --11). 

Proceeding in a similar way as in Tection it can he 
shoY^na “theft "the limits on a- ore such that, 

*— -V " 


G 


0to2Ti:, Gv^Oto 


g^G-^toG+ 5 for G > ^ 

g - Gto ? + G for . G < ^ . 

Thus the integral over R, the region in G-g space is 


I - f ^ f 

R R. R. 


where 


/ = f 


K. 


G+ C 

d G’ J" d g 
G- ? 


and 


C+G 


/' = /' f 

Rp 0 C — G 


dg . 


Rinally the integral (A-38) cai^ he written as 

TC 8 


(A-43) 


(A-44) 


(A-45) 


■ / / / (/ ■ j ) 

I’=0 9^=0 4>q_= 0 R, Rp 

- 2 

exp (- - g 2)^^(U,P,S,R,G) G^I" sin 0^ 


sin 0 


G 


dP dG dBj, de^ d4>j, dd^ 


(A-46) 



XIX 


c) Collision I^anmics and Evaluation of a (5, n ) 

I'ig. (A— 2) gives the details of a- two "body encounter 
• in the centei-of~mass system where = ~^21 “ ^ “ n 5 i i 

the unit vector along x-direction; k is the apseline vector; 
scattering angle x == Ti-2'ii defines ip , and £ is the azimuthal 
angle. 0 is the angle made by the relative velocity vector 
with x-direction. 


If and mg axe the masses of the particles whose 
velocities are respectively ^ and n the following quantities 
can be defined. 


m 


iVu 


1 


m^+mg 


mg 


and ]x 


m^ mg 

12 ~ in^jlxQg 


(A-47) 


If C* is the velocity of the particle 1 after encounter, 
then from conservation of momentum and energy v/e can write ^ 

5 ’ = 5 + 2 H 2 g cos If) k . (1-48) 

Alsu from spherical trigonometry we have the result 

k . 1 = cos 0 cos if) 4- sin 0 sin if) cos £ . (A-49) 


0 is defined by, 

a = / / (4^ - 4^) §21 ^^-50) 

b £ 

where b is the impact parameter and £ is the azimuthal angle 
in an encounter. 4i “ 4i(«£.) moment function and 

4>; = 



xiii 


He also -use the results fron ref. 18 that for a. force 
law of blie form. 


Yfe can define 


-n 2 s- 5 

®(1) _ /■‘^12s's"--'T sl'l , . 

^.0 - V-iV^ s 


i; 


(1-51) 


where 


l-l(s) = / (1- cos^x) 0.^’ 

-J- Q " ^^00 


(^^-52) 


and 


,2 s -1 






S' 3^ 


An (a) is a pure, numher depending only on 1 and s. Por 

r 1 1 

s = 2, the integrals iji ^2 said $ .^2 diverage and hence 
a cut-ofi is recLuired in evaluating these integrals. If v^.^ 
is the cut ofi value of v^, then from ref. 18 we have, 


•1 


( 5 ) = 0.422, Ap(5) = 0. ,30 


A.,(2) 

'^2 ( ^ ) 


In (1+15 i); 


In (1+-Vo-,) 


V 


.2 

o1 





( 1 - 53 ) 


(A- 53 ) 


In the present analysis Dehye radius is chosen as the cut off 
d i st an ce so t h at 

,2 


V 


o1 


9k 

g-~ 

A%e n , 



xlv 


A-III = ? 


Prom eq. (A-48) it car be written that 

= 2 M 2 g cos ijjk.i (A-54) 

Introducing eq. (/l- 54) in the eq. ( /i-50) and after using the 
result (A-49) , eq. (A- 50) can be integrated and from eqs.(A-51) 
and (a- 52) it may be virritten that 

( 1 ) 

0 = 2TC M 2 g cos 0 ^^2 * (A-55) 


a) s = 5: 

Prom eq. (,1-51), eq. (i~55) can be wnirten as 


0 = 2tc M 2 (£ 21 *^^ 'i-^(5)« 

Prom eq. (A-40) one can write that. 


1 /: 


1 


U12 


Prom eq. (A-57) , integral (a- 46) can be written as 

00 % 2% 

^ =f f f 

1-0 e^rO Pv^ 


(/ 



exp (■ 


P 


-.2 


g2) (g cos 0^-0 eV sin Bf, siD'ej 


{A-56) 


2it Ho (— ) ' ■■ Ag5) (G cos 6(j - ? )• (A-57) 




- dG-'IP d9Q d0p d4Q_ ddj, 


(..-58) 



Using the limits ±or and R 2 from (a- 4-5) and integrating 
(A-58) it reduces to, 

I = .,(5) 135 /Ri-p)V 2 , . (,,. 59 ) 

■b) 


From eq. (a-- 51)5 expression for cr from eq. (A-55) can 
he written as 


"B 

a = 2x Mg (gg^.t) (-1^)^ g“^ -1^(2). 

Hence from eq. (a-^O), 

B 

O ^ = 2tj: M 2 (~)^ ii(2) r:'^^(G cos - c ) . (—60) 


Using this -value of aid using the results ( a - 4 - 5 ), the 

integral (A-58) caii he evalua,ted and can he shown to he equal 
to 


I 


= 271 Mg (^)^ A^(2) (i_p)V2 

— exp (-?^) - Srf (?) 

_? z -I 

? = 


(A-61) 


I 


0 5 when 


0 . 



X7X 


-'■-IT (f (£.) =5 


Fiojii ref, 80 5 'the components of post collision velocity 
of one of the particles can Pe written as, 


t 


,y 




E y M2 


2(r,,,-{P oosV'^ S^-(nx-5l 


”1 

cos (£ - co^) 


2 r 2 

cos tl;+ Vg -( 

Hy- ^y) 


Sin 2^ cos (o -- 03^) 


’z 




sin 2ip cos (0 - oj^) 


whe re 




CO, 


0 , 


(n )(Ti ) 




0)^ 


( j"? y ^ ^ ^ z" _^z ^ 


[g^-(r;z-EP^ 


( 1 - 62 ) 


(A-63) 


Introducing the expressions (1-62) and (a- 63) in e q. (a- 50), 
and usina in addition to the results (1-51) and (1-52), also 
the integrals 

271; 271 

j' cos (£— w) dG = 0 5 j' cos (6— to) d£ = 


0 


0 


7t , 



x-vii 


it caji loe shown thai 


a 


2rt Og « y^(-) 




(2:v^i)a 


^-.2' 


(.1-64 ) 


From eq. (ji_40) aiid (A-8) to it can “be ^vritten 


that 


= 2 % 1\2 


( 2 M 2 - 1 )g^-U^+(U-?)^+( 1 - 2 p)r/- 2 CF cos e 


+ + (l]-Q(1-p)l Cr cos er,+2F.G 


( 1 - 65 ) 


s=5 ° 


Using the le suits from (...-51) and (A- 52) in (A-65) and 
introducing the resulting expression in the integral (A-46) , 
the integral can he evaluated and shown to he equal to 

I = (^)V2^^(5) p.V2(i.p)3/2 

■ 5(''i2“P) + 2112?^ - 2U?J . (A-66) 

s=2 : 

Using the results from (A-51) in (A-65) and suhstitutin 
the resulting exioression for in the integral (A-46) s 
the integral can he evaluated and shown to he equal to, 

I = 4it^ M 2 (—f i-,(2) (l-p)^"^^ 

^ ^12 

. -1- Urf (C) (M 2 - 2 ) 4- exp (-?^) (2 _ p) . . 

(A-67) 



xvlii 


A-V 


. 2 , 

’X" 




From a procedure similar to that in (1.--I1I), expressio 
for a can he w ritten as 

Using the results of (A-59) and (A-40) and also (.1.-8) to 
(A-IO), expression for comes out to he, 

= A-% Mg |^(U+F cos 9j,-p G- cos 0n.)(G- cos 9,^ -<;) 

+ M2(G cos ^^12^ ” "5' ^12^^ 

(A-68) 


r 2 2 

g -(G cos Gp-c ) 


1 

•3 12 J 


Using this expression for the integral (A-46) can he e’valuai 
in a similar way as done earlier for s = 5 and 2; the results 
are , 

s = 5 



.M2^^~2a^( 5)U?+ A^(5) (Mg-P)!. 


(A-69) 
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I 471^/^ pV^ (i_p)3/2 

-2pA^(2) + 2/^(2) i + (_2p..^(2) + 2M2A^ ( 2) 


- g" 412X2(2)) ^ 1 + Erf (?) A^( 2 ) 

- 12-12(2)] -I- _ 2 E^( 2 ) ^ + ( 2 p A^( 2 ) - 

+ I M2 ^2(2)) . 


2!.f2A^(2) 


(-4-70) 
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